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1. Introduction. — The first part of this paper, written by Love, contains a theoretical 
solution of the problem of rational hydrodynamics which has been named by writers 
on ballistics, (i Lagrange's problem " ; the second part, written by Pidduck, gives 
the application to ballistics. 
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In the problem* it is supposed that a given mass of gas, which is initially in a uniform 
state, is contained in a segment of a tube of uniform section. At one end the segment 
of the tube is bounded by a fixed transverse section, and at the other end the tube is 
closed by a piston of given mass, which is initially at rest and is free to move along the 
tube without resistance. It is required to find the subsequent states of the gas and 
the motion of the piston. 

Under the pressure exerted by the gas the piston begins to move, and wave-motion 
of finite amplitude is set up in the gas. The waves are plane. The theory of plane 
waves of expansion of finite amplitude has been the subject of much study, f chiefly 
in connection with the question of the initiation and maintenance of surfaces of 
discontinuity. The difficulties associated with this question do not arise in Lagrange's 
problem, because the waves that are generated are always waves of rarefaction, and 
there is no tendency to discontinuity in waves of this type. Among the results that 
have been obtained in the theory of plane waves of finite amplitude, two are specially 
important for our present purpose. The first of these is that there exist waves of the 
type known as " progressive waves/ 5 and that they are the only ones that can advance 
without discontinuity into gas at rest. They are sometimes described as " motions 
compatible with rest. "J The second important result is that the equations governing 
the propagation of waves which are not compatible with rest can be integrated. § Such 
waves will be described in the sequel as " compound waves." 

The most important writings in which Lagrange's problem is dealt with are the 
memoir of Hugoniot cited above, H. Hadamard's ' Legons sur la Propagation des 
Ondes,' Paris, 1903, and a memoir by F. Gossot and K. Liouville in ' Memorial des 
Poudres et Salpetres,' vol. 17, 1914, p. 1. 

The problem is not rendered essentially more difficult if it. is supposed that the 
segment of the tube occupied by gas is bounded by two movable pistons of given 
masses. Provision can be made for the case of a fixed end by taking the masses of 
the two pistons to be equal, for then there is never any velocity at the section midway 
between them. 

The tube will be thought of as running from left to right. When the pistons begin 
to move progressive waves set out, one from the left-hand piston with a front proceeding 
towards the right, the other from the right-hand piston with a front proceeding towards 
the left. These waves meet at the middle section, and from that section there then 
sets out a compound wave, which has an advancing front, proceeding towards the right, 
and a receding front proceeding towards the left. This wave will be described as the 

* S. D. Poisson, " Formules relatives au Mouvement du Boulet . . . extraites des Manuscrits de 
Lagrange," Paris, ' J. Ec. Pol./ eah. 21 (1832). 

f Eeference may be made to Lamb's ' Hydrodynamics,' ch. 10. 

J H. Hugoniot, Paris, < J. fie. Pol.,' cah. 57 (1887) and cah. 58 (1889). 

§ B. Riemann, 'Gottingen Abh.,' vol. 8 (1859-60); also ' Ges. math. Werke,' Leipzig, 1876, 
p. 145. 

2 B 2 
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" first middle wave/ 5 When the advancing and receding fronts of the first middle 
wave reach the pistons the original progressive waves are obliterated, reflexions take 
place at the pistons, and new compound waves are generated at the pistons and encroach 
upon the first middle wave. These waves will be described as the " first reflected 
wave from the left " (or " from the right ?? as the case may be). The reflected waves 
meet at or near the middle section, from which there then sets out a new compound 
wave called the " second middle wave." This wave again has two fronts, one advancing 
and encroaching upon the first reflected wave from the right, and the other receding 
and encroaching upon the first reflected wave from the left. The two fronts eventually 
reach the pistons, and then the second middle wave will have obliterated the first 
reflected waves, and will itself be reflected so as to give rise to new compound waves 
setting out from the pistons. These will be called the " second reflected wave from 
the left " (or " from the right " as the case may be). The motion goes on in this way 
until a piston reaches an end of the tube if the tube is of finite length. 

In what follows Articles 2-9 are devoted to giving such an account of the theory 
of plane waves of finite amplitude as seems to be necessary for the discussion of the 
problem. Although so much has been written about the subject, it appears to be 
impossible to find what is wanted in a suitable form. Articles 10, 11 contain the formulae 
relating to the two progressive waves. These are already known from the work of 
Gossot and Liouville, but it seemed to be desirable, for the sake of completeness, 
to obtain them anew. Articles 12-17 deal with the first middle wave. Sufficient 
indications of the method of determining this wave have been given by the same writers 
for the case of equal pistons. The really formidable difficulties of the problem begin 
to present themselves when an attempt is made to discuss the waves reflected from 
the moving pistons. In Articles 18-25 an approximate method of solution is found. 
It seems to be capable of giving results for the first reflected waves correct to any desired 
order of accuracy. In Articles 26-31 the second middle wave is determined. However 
far the approximation to the first reflected waves is carried, the second middle wave 
answering to them can be found by the method here given. Articles 32-40 are devoted 
to the determination of the second reflected waves. The method used for the first 
reflected waves does not give a sufficiently close approximation, and a new method is 
applied. Numerical calculation of a particular example showed that all information 
that can be of practical importance may be obtained from a solution which does not 
go beyond the determination of these waves. The results of this calculation belong 
properly to the second part of the paper. 

Theory of. Plane Waves of Finite Amplitude. 

2. General Equations. — -The motion is supposed to take place in an unlimited straight 
tube of uniform cross-section «. Let x be a co-ordinate measured along the tube, 
and specifying the position at time t of a plane of particles, which, when t = 0, is in 
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the position specified by x . At the time t — the gas is supposed to be at rest. Let 
p and p denote the undisturbed pressure and density, supposed uniform, and let p, 
p, u denote the pressure, density and velocity at time t for the particles specified by x . 
The equation of continuity is 

dx 



and the equation of motion is 



dx f 



Po> 



d 2 X dx dp 

p - * - 



dt 2 dx dx . 

On introducing u 9 which is dxjdt, these equations become 

du _ Po dp du _ 1 dp_ 

dx p 2 dt dt p dx . 

It is supposed that p is a uniform function of p, and it is convenient to introduce, 
after Riemann, a quantity a- by the defining equation 



1 



/(t) d P 



do- = - A/ [~ 
P v \dp 

and the condition that o- = when p = 0. Then the equations become 

f^ + n^- = o, — + n — = o, 

3£ 3# 3 3£ dx 

where II is a function of p defined by the equation 

The quantity II, which is of the dimensions of a velocity, may be regarded as a known 
function of o-. The value of II when p = Po is the velocity of sound waves of small 
amplitude in the undisturbed state of the gas. This will be denoted by a. The equations 
are of Lagrangian type, and x and t are the independent variables. The quantities 
p and p, like II, can be regarded as known functions of a. The value of o- when p = p 
will be denoted by o- u# 

3. Progressive Waves. — Two quantities r and s may be introduced, after Riemann, 

by the equations 

or-\-u = 2r, (j—%i = 2s, 

or 

cr = r + s, ^ = r—s. 

The equations of continuity and motion then give the two equations 

9r+n — = 0, sfc-n — = 0. 

3£ 3^a ' 3£ 3a?n 
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If s is constant, the second of these equations becomes an identity, and the first can 

be integrated in the form 

r = F (x -llt), 

where F denotes an arbitrary function. This can be proved easily, and the equation 
can be written 

X -TLt =f(cr). 

In like manner, when r is constant, the first of the two equations becomes an identity, 
and the second can be integrated in the form 

A motion with constant r or constant s is described as a " progressive wave. 5 ' A 
wave with constant s is propagated in the direction of increase of x, with velocity II, 
which depends upon the constant value of s and the local value of r . This is the velocity 
relative to the medium, not the velocity relative to the tube. Similar statements hold 
for a wave of constant r. 

4. Motion of a Junction. — When a wave is transmitted into gas at rest, or into a 
region where there is some other state of motion, there may be discontinuity in the 
values of the pressure, &c, in the two regions separated by the front of the wave. We 
consider here the case where there is no such discontinuity, but, while the pressure, 
&c, have the same values on the two sides of any plane x = const., the laws of variation 
of these quantities on the two sides of a wave-front are different. We describe such 
a moving wave-front as a a junction." Our immediate object is to determine the 
velocity of a junction relative to the medium. We shall attain this object by supposing 
that there are very slight differences between the values of any of the quantities on 
the two sides of the wave-front. 

Let w denote the velocity of the junction relative to the medium. In a very short 
time St a mass equal to p wivSt has its motion and state changed from those specified 
by u, p, p, to those specified by u + Au, p + Ap, p + Ap. The increment of 
momentum must be equal to the impulse of the difference of pressure, and therefore 

we have the equation 

Pqww St Au = oo Ap St. 

Further, the work done during the interval St by the external pressures on the ends 
of this element of mass must be equal to the sum of the increments of the kinetic and 
intrinsic energies of the element. Now the changes of state being adiabatic and very 
slight, the increment of the intrinsic energy per unit of mass may be put equal to 

~pA(l/ P \ 
and therefore we have the equation 

w(p + Ap)(u + Au)St—a>puSt = %p oowSt {{u + Auf-ic 2 } —p wwStp A(l/p). 
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The two equations containing Ap and A^ give 

Ap = p w Au 
and 

p Au + it Ap = p Q wu Au + p w {pjp 2 ) Ap. 

The terms uAp and p wuAu in the second of these equations cancel, and then, by 
eliminating Au between the two equations, we find 

p Ap 

Since there is no actual discontinuity and p is a uniform function of p, we may replace 
ApjAp by dp J dp, and thus obtain the equation 

p dp 

which shows that the velocity of the junction relative to the medium is that which 
was previously denoted by II. 

If motion is set up in one part of the gas, and advances into previously undisturbed 
gas, the value of p at the junction is p , and therefore the velocity of the front of the 
wave, relative to the medium or to the tube, is that which has been denoted by a. 

5. Nature of the Motion in a Compound Wave. — Important results can be obtained by 
regarding x and t as functions of r and s. On interchanging the dependent and inde- 
pendent variables in the equations 

3r , t-!- dr A ds -^ ds A 

— +n — =. o, — -II — = o, 

at ox at ox 

we obtain the equations 

cXq tt ct __ n ox . tt06 _ n 

3 XX _ V/, ^ I XX ^ \J . 

6* as or or 

Now the differentials of x and t are always connected with those of r and s by the 
formulae 

7\rp GOG dt ut 

dx = -t-^- dr + -t~ ds, dt — — dr + — ds. 

dr as or os 

Hence the places in the medium, and the times, at which any particular value of r is 
found, vary according to the formulae 

dx = |k ds = n | d s> dt = f 8 <h = jj[> 

and thus it appears that any value of r is transmitted through the medium, in the 
direction of increase of x, with the velocity II. In like manner it can be shown that 
any value of s is transmitted in the opposite direction with the same local velocity. 
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We have seen that the velocity of a junction relative to the medium is the value of 
II at the junction, and it follows that the value of r remains constant along a junction 
which travels in the direction of increase of x. If the junction travels in the opposite 
direction, the value of s at the junction remains constant. 

The motion consequent upon any initial conditions consists in the transfer of the 
existing values off and s through the medium with the variable velocity already described. 
New values of r and s can be generated at boundaries and transferred through the 
medium. 

6, General Analysis of Compound Waves. — -When the dependent and independent 
variables are interchanged in the equations 



dt ' dXn dt ' dxc 



there result the equations 

ou ocr da- ou 

The first of these shows that there exists a function Z of a- and u which has the properties 
expressed by the equations 

— , t = — - 5 

do- OU 

and then the second shows that Z satisfies the differential equation 

— ( : n —\ — tt ^ 2 ^ — o 

If Z can be found in accordance with this equation, the values of x and t answering 
to any simultaneous value of a- and u can be deduced. 

There is a relation between Z and x 9 which can be obtained very simply by introducing 
for a moment a quantitv m by the equation 

& = Po/p> 

for then we have 

__ do- 
ll 

and it follows that we have at once 



&dp = dm, 

P 



v£a , ■ v/j 

dm ou 
and 

— Qua. ' — ^ 

OJU(\ ot> 
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These are the relations of duality familiar in discussions of partial differential equations,* 
and we may put 

ZOOb , . OJb 
— ry* _1 f _____ -■■ f"Y> 

OX dt 

Actually Z could differ from the right-hand member of this equation by a constant, 
but as such a constant would be irrelevant, the above will be taken as the relation 
between Z and a?.. 

The equation satisfied by Z can be written in either of the forms 

ccr 2 \II dcr J dor du 2 



or 



3 2 Z / 1 dll \fdZ dZ\ __ 



3r9_r \2IT dar/\dr ds 

7. Relation between Pressure and Density.— -The analysis of the problem is not rendered 
more difficult if the adiabatic relation between pressure p and volume i? is taken in 
the form p (v — 6) Y = const, instead of the more ordinary form jw y = const., and 
the former is more suitable for the applications which we have in view. We shall 
accordingly take the relation between pressure and density to be 






1 IV (I IV 

\p pi \p pi 

where /3 and y are constants. Then the following results can be obtained without 
difficulty : — 

- JL_ f Pov(P-Po) y ( fa, P Y Y " i)/8 

y— 11 fipQ J \ p fi — p, 

v = 2 [ Poy(iQ-Po) l^ 

y— 11 /3/> 



(7 = 



a 



Po(P-Po)} ' 



P=Po(<rf<r Q )** +1 > 

Il = a (<r/<r )* 

where 2w has been written for (y + I)/(y — I). 
The equation for Z can now be written 

3 2 Z , 2n 3Z 3 2 Z 



or 



1 _____ __ = o 

<r" a- da ou 



H - ^ r tt~ = U. 



dr ds r + s \3r * 3 s / 

* The reduction of the equations governing the propagation of plane waves of finite amplitude to a 
single partial differential equation of the second order was effected by Riemann, who worked with " Eulerian " 
equations. The use of the principle of duality to connect Z and x was noted by Hadamard. 
VOL, COXXII. A. 2 
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8. Integration of the Equation in Special Cases.— When n is a positive integer, the 
equation can be integrated. We write for a moment D for d/du, and observe that, if 
D were a constant, the equation 



3 2 Z 2n 3Z T ^ 2y 

oar a ccr 



would be a form of Riccati's equation, and could be integrated in the form 

1 aV-M^A + e-^B) 



.0- C(T 



O" 



where A and B are independent of a-. Treating them as functions of u, we obtain the 
general primitive of the equation for Z in the form 



Z = 



1 IV" 1 \ F(<r + u)+f{cr-u) \ 
,<r do-/ I cr j 



9. More General Integration. — Interpreting the variables r and s as the co-ordinates 
of a point in a plane, Riemann showed how to integrate the equation for Z when the 
values of this function and its first differential coefficients are given along an arc of a 
curve in the plane. If V satisfies the " adjoint " equation 



3 2 v / a a 

—n{ ^~ + 



the integral 



3r ds 



3_ fy / 3Z nZ 
or\ \ds r + sj 



V N 



3r * 36V \r + s, 



os { \ or 



= o, 



nV\\ 
r-\-s/j_ 



dr ds, 



taken over any area in the plane, is equal to 



V<!;£i- + 



n 



'dZ dZ" 



drds r + s\dr ds / 



Z 



3 2 Y /3 ( 3W V \V 
3r3s \3r 3s/ \r-fs/J_ 



dr c?s. 



and therefore vanishes. It follows that the line-integral 



J \ ds r + 5/ 



'3v___^v; 

3r r + s/ 



Adr 



taken round the boundary of the area vanishes. 

Let the values of Z and its first differential coefficients be given along an arc AC of a 
curve, and let P be a point which is not on the arc. Through P let lines PA and PC 
be drawn parallel to the axes of r and s 9 and let the area of integration be that bounded 
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by the arc AC and the lines CP and PA. 
integral is 



The contribution of CP to the line- 



cp \os r + s/ 



which may be written 



[VZ] P -[VZ] 



c 



Jla 



'av 



- ) ds. 



CP 



3,9 r + sj 



The contribution of PA to the line-integral is 

pa \ or r + sj 




Y 



Fig. 1. 



Now we can find a function to satisfy the equation for V, to make V = 1 at P, and 
so that, along CP, where r has the same value as at P, dY/Bs = riVftr + s) 9 and along 
PA, where 6* has the same value as at P, SY/dr = riVftr + s). Then the value of 
Z at P is 



[vz] 



ac I \os r + s, 



dr r + s) 



)dr. 



The required function V can be shown, after Riemaot, to be given by the equation 



V = ^)"F(n,l-n,l,a 



where F is the symbol for the hypergeometric series, 



i 



(r-V) (s— s') 

(r + s)(r f + s f )' 



and /, s' are the co-ordinates of P. 

It may be observed that if n is an integer the series terminates, and V like Z is 
expressible in a finite form. 

It will be useful hereafter to note the formulae 

9V nV = (r + sY-tjs + ^is- s') d v( , 

dr r + s (r' + s') n+1 d£ [7h n > L > V> 

9V _ nV_ = _ (r + s)"- 2 (r + $')(r-r') d v( ^ 

35 r + s " (S + s')* +1 <# ^ 

The Progressive Waves in Lagrange's Problem. 

10. TAe Progressive Wave from the Left— Let the positive sense of the axis of a; be 
from left to right, and let the initial positions of the two pistons be given by x = 
and ;r == c, where c is positive. We shall denote the mass of the piston at x = by 
M, and that of the other piston by m. 

2 o 2 
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The progressive wave generated at the piston M is determined by the equation of 
motion of this piston. This equation is 

M -^— = —Mp, 

ct • 
and it must hold at x = for all positive values of t. It may he written 

M^ = -«p (-J , . 



and, since in the progressive wave s is constant and equal to |-o- 0? or <r — w = <r , it 



gives 






wp 

from which, since o- = o-,» when £ = 0, we have at # ft = 



JxL(7n I / <T \ 1 L 

1 i * 



2fla)p [\<r } " ) 

Put for brevity 

H = Mo- a/2nwjp , 

then we have the values of <r and £ at # = connected by the equation 



a [\<r 



Now in the progressive wave we have 

x Q -IH =/(«■), 

where the function /is to be found from the condition that at x () = the above relation 
holds between cr and t. Hence we find 



/(,)=-Hji-^n- 



l \<T / J 

and the progressive wave formula can be written 

at + B. /V 2 * 



o 



a? -f H \cr 

In the motion described by these formulae any plane of particles, specified by a value 
of x in the interval \ c > x > 0, remains at rest until t = a5 /a, and then moves with a 
velocity u, which is equal to cr — n> Therefore the value of a; answering to these 
particles at any subsequent time is given by the equation 
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X ' — Xa "r 



u dt = a? + (oo — <0 2 ^ — — - ~~-7 do-, 



,2n-fl 



or 



«/? — - .^/q "T" \ Xq ~j~ JUL ) 

a 



2n 



2^i — l 






2n-l 



I 



<*0 
,0* 



2w 



This equation holds so long as the plane of particles is in the region occupied by the 
progressive wave. In particular, the displacement of the piston M is given by the 
equation 



x 



2n 



— 7 (o- — "■) — + ( or-cr* t 

— 1 a \2n — l / 



in which 



t = 



H 



cr. 



-° 



2» 



a Ivr 



V. \«-» / 



The corresponding values of Z are found from the formula 



in which 



to be 



dx 


„ Pn 


dx Q 


p 


HJ 


2n 



(jX 

Z = x - Yut—x, 

OJur\ 



/3—p p 



cr 



2n~l 



Po /^ — p \<r c 



o*. 



\2w-l" 



a l2n— 1 



(7a — or- 



o 



2w,— l \cr,/ 



11. The Progressive Wave from the Bight. — The equation of motion of the piston 

m is 



and we put 



du 
m — = «>p, 



/^ = mcr a/2nwp . 



Since in the progressive wave r is constant and equal to J<r , or o- + m 
of cr and £ at # = c are found to be connected by the equation 



2n 



t 



h I / (To j __ -j 

a l\or 



,1 



O"0: 



the values 



and then the progressive wave formula is found to be 



at + & 



CTr 



2n 



c-\-h—x n \<t 



The value of x for any plane of particles specified by a value of x in the interval 
Jc < # < c, and for any time later than that given by t = (c — x )ja, is found to be 
given by the equation 
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a 



rr \ 2n 
^0 ) 



2n 



<Tr 



2»-l 



I 



or / 






2n-l W 



which holds so long as the plane of particles is in the region occupied by the progressive 
wave. In particular, the displacement of the piston m is given by the equation 



2n ( \ h I 2n 

2n — 1 a \2n~-l 



in which 



£ 



M, 



fl9 ) 



2» 



flf — <T U. 



U- 



The formula for Z is found to be 



z=, c+/ * JhT 1 -! 



{2n-l)a 







\ o- J 



h 



a 



(cr — <r). 



The First Middle Wave. 

12. Conditions satisfied at the Receding Front. — In the progressive wave from the 
left s is constant and r variable. The greatest value of r, which is the undisturbed value 
■Joro, travels at the front of the wave, and continually diminishing values of r, generated 
at the piston M, travel after it. Similar statements, with r and s interchanged, hold 
for the progressive wave from the right. The fronts of both waves travel along the 
tube with velocity a. When they reach the middle section, a compound wave begins 
to be generated there, and transmitted in both directions, encroaching upon the original 
progressive waves. This wave has a receding front, along which s is constant, travelling 
towards the left, and an advancing front, along which r is constant, travelling towards 
the right. The constant values of r and s at the two fronts are equal, and each of them 

IS -^(Tq. 

At the receding front the variations of x and t are connected by the equation 



dx + TL dt = 0, 

while the values of x , t and a- are connected by the progressive wave formula, which 
can be written 



tf -m+H fi- 



ll 

a 



0, 



so that the variations of x Q , t and II are connected by the equation 

dxo-{ndt + tdm- — dIl = 0. 

a 
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On elimination of dx there results tlie equation 

\ a 
which can be integrated in the form 

HY 

M ) II = const. 

To determine the constant there is the condition that when x — fc and t = %c/a, the 
value of II is a. Hence at the receding front we have 

or 

(,z' + H) (at + H) = (|c + H) 2 . 

13. Conditions satisfied at the Advancing Front. — At the advancing front we have 
in like manner 

1 + — J = 0, 

leading to 

(at + hfll = (%c + h) 2 a 
and 

(c + h—x )(at + h) = (ic + A) 2 . 

14. Conditions determining the First Middle Wave. — It will now be convenient to 
restrict the value of n to be an integer. This happens when y has one of the values 
3, 5/3, 7/5, 9/7, 11/9, . . . With a view to applications, in which the value of y 
is 1 • 2 nearly, we shall take the value 1 1 /9 for y, or 5 for n. Then in any compound 
wave Z has the form 

or 

105^ 9 F(^+tO-105^ 8 F (1) (<r+tO + ^^ 

where F (1) , F 00 , and so on stand for the first, second, &c, differential coefficients of 
the function F with respect to its argument. We have to determine the unknown 
functions from the values of Z at the advancing and receding fronts. 
At the advancing front, where r = |o- and a- = |o- + 8, we have 
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and at the receding front, where s = |o- and o- = |-o- -f • r, we have 



rr 



ITT ^0 



9 



<T, 



o 



a 



9 1 TT 
i<r + r/ J a 






15. Determination of the First Middle Wave,- — To determine Z from these condition 
we may have recourse to Riemann's method, taking the curve AC to consist of 

segments of two lines AB and BC, which meet at 
the point B, where r = s = |o- , and are parallel 
5 respectively to the axes of s and r. 

We have then 




A 



tt /oZj 5Z 
cp \ds r + s, 



CM +1 Z ( — — 

pa \3r r + s. 



r 



Fig. 2. 
This equation is 

[VZ] P 



/a 7 p;7 
ab \ 3s r + s, 



r ,/av 5 V 



BC' \ 3/ 



T + S, 



c?r 



d!r = 0, 



[VZ] 



c 









V 5V 



+ [VZ] B -[VZ], 



jcp \ 3 s r + Sj 

rr / 9V 5 V 

^ 

. ab \ 8s r + s, 



C?S + 



ds + 



PA \ O'?' 



''3Y 5V 



•r + 67 



) c?r 



^" 



BC 



Z; U V J V 
t TT ~ 

or r + ^, 



c?r = 0, 



or 



Z (/, J) = [VZ] A - [ VZ] B + [VZ] C + 



• „m 5V\, 

Z -~ ds 

ab V os r + s 



■ 7 /dY 5V\ , 
bc \ or r + sj 



where /, s' are the co-ordinates of P. Also at A we have 



r 



<T{\* S 



2 U 0? 



«', £ = o, V 






z 



| (c + h) 



(Tr 



a 



or. 



9 







so that 



[VZ] A = -|( C +A.) 



er. 



o J 



CTr 



9 



At B we have 

so that 

At we have 

/ — * / . o - — > *> u. 



2°0? 



so that 



C"VZ] C 



a l(r'+V) 5 (|<x + /) 4 



r = <s 



2 (T 0? 






Zf = (), 



vz 



B 



0, 



£ = 0, V 






r -t-s ' 



£k == 



xl — - 



J) 



a 



,2^0 + ^ 



l\-^(fcoS): 



Nko-^Ofe-rT' 



a 



r + s 









] XT ^0 
tr JLx 



<r, 



9 



o 



9 



>l(r'+/) 5 (K + ^) < 



^rt + r'^ 5 



zrlK-f^te' 5 



IT 
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Along AB we have 



a 



arc 



\9 



<Ta 4- S, 



>2°o-r 






where 



8V 5V 

3 s ?* + s 



(l^ + s^^^r^^^ (20-180f +.420f a -280f 8 ). 



f = - 



so that 



(W-r')(s-s') 



5 






■ids = 



i"a 



*s' 



(c + A) 







0", 



9 







9a 



,2-tTo + ^. 



1 









x 



20 (K + ^) 3 + 180 fez^ ^~^ (^o + s) 2 



r + 5 



Along BC we have 



+ 420 (2^)-^) 2 (^-^) 2 (^o + ^) { oo (k-^) 3 (^-^) 

(r f + s') 2 (r / + .9 / ) 3 



3" 



cfc. 



7 — — XTT 



0"f 



a 



i^o + r, 



9 1 TT 

J a 



3X _ AY_ _ (K + rQ&ro-sQfero + r) 8 



3r 



T -f ,9 



where 



(r' + s') 6 



(20~180<; + 420£ 3 -280£ 3 ) : 



& 



so that 

rz /av 



(r' + s'H&xb + r) 



BC 



O 



"5 V \ , 

- a? 1 



dr r + ,s7 



2'0-0 



9a (Af<x + r/ j a 



(K+^)(K-^) 

(r' + s') 6 



20(|-(7 +r) 3 + 180 



r' + ,s' 



+ 420 fezV)l(r-^(^o + r) g() 

(r + s ) 



(V~fV) s 



?r, 



a: 



There is no difficulty in" the integrations. After they have been performed, we may 
suppress the accents on r' and s\ and so obtain the formula for Z in the first middle 
wave expressed as a function of r and s. We find, for example, as the coefficient of 
(c + h)cr j9a the expression 



\9 



/l 



K r + s, 



fro ± g V _ (i<ro + s) (ko-^) 



; 



r +.5 



(r + s) 6 



iW-5(K + *)' 



(If^^^ 



(r + s) 

(r + s) 8 



{l89<r u 4 -420<r 3 (K + «) +270O-0 8 (|^ + .9) 2 -35 (io-o + s) 4 } 



^ ±^ ^=^{126^-420^ (K + s) + 540a 2 (|<r + s) 2 -315«x (£*„ + *)< 



+ 70(£<r + s) 4 }. 
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In this expression we put 

th>-r = ( : k<r + s)--(r + $) = i(<r + <r-u)-<r 9 

and find that the expression is the same as 

\<r 3o"/ [ o- J 

where 

r i 

<Pi {<r-u) = ro-5 °"o°- ^(o-o + o— "^) 5 {I26<r 4 -210flr 3 ((r + a—^) + 135<r 2 (<r + o— w) 3 

The remaining terms in the expression for Z may be treated in the same way, and. we 
obtain finally, as the expression for Z in the first middle wave, 

9a \ar car/ \ <r ) a \ar car) \ or j 

Hop /l 3VI^i(q-m)l H/i ayj^^ + w) 



9a \<r 3cr/ I a- J a \ar da-. 



ar 



in which the expression denoted, by c/>i has been written down, and \}s x is given by the 
formula 

ir x (<r-u) = - Y^2"9 (o-o + cr-^) 6 ((To-flr + lO 5 - 

It may be observed that the differential coefficient of the function <p x . is given by 
the equation 

o 

Although the actual calculation of Z is rather long, it is comparatively easy to verify 
that the form obtained satisfies the conditions by which. Z was determined, 

16. Transformation of the Formula,— The form taken by Z in the first middle wave 

is 

\<T 3fT/ ( (T J 

where 

H FT 

9a a 

so that <J>i and ^ are rational integral functions of the tenth degree. Now it is impor- 
tant to observe that when 1 •< v < 8, 

idViicr+uYi , lV /i ayj(<r-w)'i 
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while f or v = 9 we liave 



I 3 V \k±^\ = 2 (2 . 4 . 6 . 8) - (i # Y \lz=2£ t , 

,<T OCT J I <T J \<T COrJ I <T 



and for v = 10 we have 



10 



■1 ayj(^V!U 2tt ( 2 .4.6.8. 10) + (I AY ji^l 

,a- coy L <r J v \<r do-/ I o- 

It follows that the expression for Z can be written either in the form 

Z = (^J{iM} + K 1+ L A 
or in the form 

\ar Car! I o~ J 

where 

Kj = 2 . (.2 . 4 . 6 . 8) x the coefficient of {a-—uf in < \r l (ar—u), 

Lj = — 2 . (2 . 4 . 6 . 8 . 10) x the coefficient of (a-—u) 10 in ^ (o— u) 9 

h i = 2 .(2 .4 .6 .8) x the coefficient of (o~ + ^) 9 in $ x (<r + u), 

l x = 2 .(2 .4 .6 . 8 . 10) x the coefficient of (a- + u) 10 in 4^ (o--fu), 

and Qi and j x are certain rational integral functions of the 9th degree. The explicit 
expressions are 

Kj = ~l (c + h) (a-Ja), L x = — h/a, h i = |Ho- /a, l x = —Ufa, 

+ ^^2^ { 3 1 5,r 8 (o- + w) - 420,r 6 (<r + «) 3 + 37 8<r 4 (cr + m) 6 

-180<r 2 (<7 + w) 7 + 35(<r + tt) 9 }, 

and g x (o— w) is obtained from Q, (<r + u) by writing — (a-—u) for (cr+w). 

17. Incidence of the First Middle Wave on the Pistons .—The values of all the quantities 
at the piston M, at the instant when the first middle wave reaches it, are to be found 
from the formulae, connected with the receding front of the wave, by putting x = 0. 
We see that the receding front teaches the piston M at the time T u where 

T = (ic + H) 2 H , 

that the corresponding value of a- is E x , where 

v - ( H V/5 
^ - "-o \i c + H , 

2 D 2 
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that the corresponding values of f , s, u are 'R u S l5 U u where 

and that the corresponding value of Z is Z l5 where 



7 — ±± Jiiin- __y 



In like manner we see that the advancing front reaches the piston m at the time 
t u where 

t ... (h+W h 

and that the corresponding values of <r, r 9 s 9 ^ 5 Z, are <r u r u s u u l9 z u where 



The First Reflected Waves. 

18. Conditions determining the First Reflected Wave from the Left. — After the instant 
t --== TV the formulae belonging to the progressive wave from the left cease to hold in 
the neighbourhood of x — 0, and a new compound wave, the first reflected wave from 
the left, is generated there and encroaches upon the first middle wave. The junction 
is characterized by the value R x of r. The conditions determining the reflected wave 
are the condition which holds at the junction, where r = ll u and the condition which 
holds at the piston, where x — 0. It is further necessary that x should vanish when 
r = R L and s = Si. 

The condition which holds at the junction is that the value of Z, calculated from 
the formula? belonging to the reflected wave, should be equal to that calculated from 
the formulae belonging to the first middle wave when r = R^ The condition which 
holds at the piston is the equation of motion of the piston, viz., that at x = 

The condition that x {) should vanish when r = ll x and s = S x is the condition that 
dZ/do- should vanish for these values of r and s. 

To express the condition which holds at the piston in terms of Z we substitute 
jj (<rf<r ) n for p, and 

dx {} f/dxp St ^ Sxq dt 
dor I \3cr du Bu dcr 



ON LAGRANGE'S BALLISTIC PROBLEM. 187 

for dujdt. Then we liave 

M dx () __ ____ /^_\ n /8^o 3£ _ ^Ep^ \. 

cop Oct \crj ' \d<r dll du do*/ 

Again we substitute — II 3Z/3o- for x and dZ/du for t, and put 3Z/3<r = 0, obtaining 

the equation 

8 2 Z 3 2 Z / 3 2 Z \ 2 = 1 Ollo-o 10 3 2 Z t 

8cr 2 3i// \dardu) a<r U dcr 2 

The condition which holds at the junction is that 

for all values of a- and w for which & + u ~ 2R X . 

19. Determination of the First Reflected Wave from the Left. — These conditions can be 
satisfied by assuming for Z the form 



Z = (-7T- ) i-^ L \ + K 1 + L 1 u, 

\(T OCT/ I Or J 



expanding the unknown function Fi in the series 

Fi(<r+^) = A + A 1 (a-+?^-2R 1 )-fA 2 ((7+?^-2R 1 ) 2 4- .-.., 

and finding the coefficients of this series. 

The condition which holds at the junction determines the coefficients A , A 1? ..., A 4 . 
The condition that 3Z/3cr vanishes when r = Ri and s = Si determines the coefficient A 5 . 
The remaining coefficients are to be determined by the condition which holds at the 
piston. 

We have 

z = 105 Iik+A _ 105 mz+A +45 ^±^1 

(T (T (7 



-io ^(^) + ^ w ^) +Ki+L 

from which we find 



6 5 ' --i ' ~~ lU > 



A u = Q 1 (2E 1 ), A, = Q 1 (1) (2E 1 ) 5 2! A a - Q.^E,), 3! A 3 = Q 1 (8) (2K 1 ), 4! A 4 = Q 1 (4) (2K 1 ). 
We have also 

8Z = _ 945 j^«) +945 F^(o ;+M ) _ 420 j^k±tf) 

Oar <r w cr <x s 

0" cr cr° 
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and the condition that this vanishes when r = R x and s = Si gives 

-945^+945^-420^^ , 

^i ^1 ^1 ^i ^i ^i 

thus determining the coefficient A 5 . It is seen easily that 5 ! A 5 = Q a (5) (2K!). 
Now when dZ/da- = the differential equation for Z shows that 

c Za d Zi 

dcr 2 du 2 

Also we have in general 

8 8 Z __ 1Q5 jVlktiQ 105 F 1 (3) (q-+^) i5 F 1 w (.+ w ) 10 F 1 fe) (^+^) , F^Qr*^ 

(TO or o - cr cr cr 

•K = - i ,45 F '"'r, + " > +a 45 F ''"'- + " > -420 F -'"'r" ) 

9<r3tt a- 10 cr 9 cr 8 



a cr cr' 

and therefore when 3Z/3o- = we have the equation 

~{945F 1 (1) ~945o^ 2) + 420^ 

+ 10H(oV7^){l05F/ 2) -105o^ - 

as well as the equation 

945F 1 -945o~F 1 (1) + 420o-^ = 0. 

The equation expressing the condition which holds at the piston is linear in F/% 
and therefore can be solved for Fj (0) without ambiguity. As it holds for r = Rj and 
s = Si, it determines the coefficient A 6 . The equation in question holds for all values 
of cr and u for which the equation expressing the vanishing of dZ/dcr holds, and it can 
therefore be differentiated totally with respect to o\ u being treated as a function of a- 
in accordance with the equation 9Z/8er = 0. This process yields an equation which 
determines the coefficient A 7 without ambiguity. A second differentiation yields an 
equation from which the value of the coefficient A 8 may be found. By proceeding in 
this way we may obtain as many of the coefficients A as we wish. 

This method of determining the coefficients A.., A 7 , ... is not very well adapted to 
numerical computation, and other methods will be explained presently. 

20. Determination of the First Reflected Wave from the Right, — The junction of the 
reflected wave and the first middle wave is characterised by the value ^ of s. The 
conditions determining the reflected wave are the condition which holds at the junction, 
where s = s u and the condition which holds at the piston, where x = c. Further, 
x must be equal to c when r = r x and s = s u 
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At the junction, where s-= s u the value of Z must be the same whether it is found 
from the formulae belonging to the reflected wave or from those belonging to the first 
middle wave. At x = c the equation of motion of the piston, viz., the equation 

cu 
m — = w p 

at 

must hold. The equation — II 3Z/3<r = c must hold at r = r x and s = s^ 
To express these conditions it is convenient to write 

\<r 3<r/ \945<W 

then Z ; satisfies the same differential equation as Z, and dZf/da- vanishes when x = c. 
The condition which holds at the junction is that 



\<7 00-/ [ cr ) 



for all values of or and w for which a- — u = 25 1# 
The condition which holds at the piston is that 

3cr 2 3?/ \3o-oW ao- 11 3a- 2 

when c7J jd<r = 0. 

These conditions can be satisfied by assuming for Z the form 



Z = (- 



Al 3 



\4 



3, 



cr dcr 



r ^ 10 



£&'■<'->. 



4- hi + ^, 



expanding the unknown function f x in the series 



f\ = ^0 + ^.1 ( cr -~^-~2,9 1 ) +a 2 (cr~?/-~25 1 ) 2 + .... 



and finding the coefficients of this series. 
For the coefficients a , a l9 . . . , a 4 we find 



CO", llj 



a ° + ^t; = ^ (2*0, ^ - g, (1) (2/0, 2 ! a, - q® (28,1 3 ! a, - g® (2*), 4 ! a A = ?« (2^). 



945a 
The coefficient a 6 is given by the equation 

_945^ + 945^-420 2 -~* + ^ 0. 

oy u ay ay oy . cry cr/ 

The remaining coefficients can be determined from the condition which holds at the 
piston in the same way as the corresponding coefficients in the formula belonging to 
the first reflected wave from the left could be determined. 
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21. Relation between Pressure and Velocity at a Piston. — The equation in terms of 
cr and u, which holds for the first reflected wave from the left" at x = 0, is the relation 
between the pressure on the piston M and its velocity during the time that the wave 
is being generated. It may also be interpreted as the equation of a certain locus in the 
plane of r and s. This locus passes through the point (R l9 S^,. and we may take its 
equation to be of the form 

r-R, = B 1 (s-& l ) + (B a /2 1 ){s-$ 1 ) 3 + {Bj2 1 a )(s-& i y + ... • 

Now if the coefficients B were known, we could determine x , as a function of r and s, 
from the known value, zero, of the function along the locus and the values of its differen- 
tial coefficients along the same curve. These differential coefficients also are known 
along the locus. To prove this and obtain formulae for these differential coefficients, 
we write X for x and observe that the equations of Article 6 show that X satisfies 
the differential equation 

i/j.3X\ "JL/I ra\ 

duXIIduJ 3cr\n 3o-A 
which can be written either in the form 



a 2 x 10 ax a 2 x 



2 ^ . ^ 9. v ? 



or in the form 



dor* or da d'lff 



0. 



a 2 



b 6A. , oA 



— - + -7T- = 0. 



3r 3s r + 5\3r 3a 



Further at # rt = we have 



£H — — ^P — M Pi)f (T \ n _ a<T o /<r Xl1 



and 



3^ M M V / 10.HW 

3a3 rr 9*o /cr\ 10 3x 



n ^> . a 

\ /T 

0' 



a?i 8cr " dcr \cr n / 3cr 



rmrutimrnritiH i n ams mmmsa mm m ^mum* —«-» -~».....» -..■„.».— a ~ 1|T —_ . — . — ™~ rm rutv-fmmmmmmm — i 



3cr 3-w 3n 3<x [dcr; \du) \3<r/ [' \d/u/ , 

where dcr/ab is to be found from the equation connecting r and s. Thus we have along 
this locus 

1.0H — 

3a? __ 10H a^ dn 



The equation for X is similar in form to that for Z, and may be solved by Riemann's 
method. When this is done the coefficients B in the equation of the locus may be 
determined by identifying the values of X at r = R x with those given for x by the 
formulae for the first middle wave. 



ON LAGRANGE'S BALLISTIC PROBLEM. 



191 



22. Integration of the Equation for X. — We write the equation 



a 2 x 



5 / 8X 3X\ _ A 



87-85 r + 5\3r ds 



and consider also a function Y which satisfies the adjoint equation 



d ' Y + 5^ + d 



dr ds 



Y 



\dr ds/ \r + $/ 



)=0. 



Then the integral 



j j 



a iwax 



87' 



Y 



8s 



_5X 

r + s. 



ds 






aY- 5Y\ 

+ • ) 

r + sj 



dr ds 



taken over any area in the plane of (r, s) 
vanishes, and therefore the integral 

y/p-^Wx^ + iX-W 

\ 85 r + s/ \or 7- + ,§/ 




taken round the boundary of the area also 
vanishes. 

We take the area of integration to be 
bounded by an arc of the locus along which 
X = 0, and two lines parallel to the axes of s and r and meeting at the point P, where 
r = r f and s = s ; . Let these be the lines PA and PC. Then we have 



Fig. 3. 



[YX] A -[YX1 



pa \ os r + 5/ 



AC 



Y i 



\ ds 7- -hsf \cr r + 67 

/8Y 5Y 
X ( 1 ) dr 



or, since X = on the arc AC. 



CP 



dr r + s, 



0, 



[YX] P = 



Y -7—- c?5 — 
ao 8s 



pa \os r + s, 



CP 



X ( — — l dr. 

\ or r + s; 



We choose Y so that, at P, Y = 1, along PA, where r = r f 9 dYjds = — 5.Y/(r + s), 
and along CP, where s = s', 8Y/8r = — 5Y/(r + s). Then the value of X at P is 



Y^rls. 

AC OS 



The requisite form of Y is 



Y = 



r f + s /x 5 



F(6, -5, 1, f), 
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where 



and we have 



(r' + s') {r + s) 



X (r\ s') 



V /^ __ Q%:\ 1 (drr — fliA 

AC \C(T dll I 



= J 



Y 

AC or 



10H du 



'dcrV 2 

du) 



Ida 
\du 



1 )du 



= -5H 



Y 



g5w. 



acct 



23. Determination of the Coefficients B. — The integral (Y/cr) ^ may be evaluated 

Jac 

approximately by assuming, as in Article 21, that the equation of the locus, of which 
AC is an arc, is of the form 

where § stands for (s — SJ/Sj. Then along AC we have 



u 



Y 



S 1 {(B 1 -1)^+B#+B 3 ^+... 

S x {(B 1 -1)+2B 3 5+3B 8 ^+...} <M, 

(r ; + s') 5 



(7 



i +30 (7r° ( r g/) +2io. ^~r /)2 , ( rf + . . . 

(r-f s')cr (r +s ) 0" 



Also any inverse power of a- can be expanded in powers of § by means of the equations 



cr 



which give 



o- = 2 L {I + (B l +l)S+B 2 P+B 8 P +...}, 

- 2i l 1 ~' f ~sr + -¥r-l-srj— -^ 



— /c 



l- ff (B 1 + l)5-^B 2 -%tl)(B 1+ i)4^ 



2! 



/cE ;3 - g ( g +l)(B 1 + l)B 2 + g(g+1 3 ) ! (,c+2) (B 1 + l)4^ 
- {*B 4 - ^1) B/-, (*+ 1) (B 1+ 1) B 3 + ' c( ' c+1 2 ) ! (g+2) (B 1+ 1) 2 B, 



.<c(/c+1)(ic + 2)(k + 3) 

4! 



(B x + 1)4 <**-... 



To obtain the expression for X on r = Rj we have to put B x for r', so that 



r-r' = S 1 (B 1 «5-fB./+B 3 <f +...), 
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and for s— s' we have to put ^(S— §'), where $' stands for (s' — S a )/2 X . If the expan- 
sions are carried as far as the fourth order, the result is that, to the fifth order in $', 



(Y/o-)du = (B i -l)§ f +{B 2 -{SB 1 -2)(B 1 -l)}S' 2 

+ {B 3 -3B 2 (2:^-1) + (7B 1 2 -6B 1 +2) (Bi-1)} S' 3 
+ {B 4 -2B 3 (3B 1 -l)-3B/+jB 2 (42B 1 a -33B 1 + 5) 

-i(28B 1 s -2lB 1 2 + 9B 1 -2)(B 1 -l)}^ 
+ {B 6 -lB 4 (30B 1 -7)-6B 2 B 3 + T \>B 3 (210B/-108B 1 + 7) 

+iB 2 2 (l05B 1 -27)- 1 VB 2 (560B 1 3 -462B 1 2 + 8lB 1 -7) 
+ |(l26B 1 4 -56B 1 3 +2lB 1 2 -6B 1 + l)(B 1 -l)}<S /5 . 

Now at any point (B l5 s') on r = B x the formula for the first middle wave gives 



0Cc\ — "— " Ob 






10 



3 



'0 



cr, 



0' 



a 



cr 



1 avrQ^cr+uy 



a 6 or J 



cr 



+ K 1 +L 1 t& 



or 



a { -945Q : (2Ej) + 9450-QW (2^) -420<rU (2) (2E0 + 105o- 3 Q 1 (3) (2BJ 

- 1 5<r*Q» (2R,) + o-'Q® (2R0}, 



10 



in which we have to put a = S 1 (l+5 / ). Then, since a? vanishes with <5 ; , we have 
without any approximation 



•A/ 



= _ JL [{945Z 1 Q 1 (1) (2E0 -8402/Q® (2E X ) +315Z 1 3 Q 1 <3) (2E,) -602/Q^ (2E X ) 

+ 52 1 6 Q 1 < 5 >(2E 1 )K 
+ { - 4202 1 2 Q 1 (2) (2E 1 ) + 3 1 SS^Q^Ri) - 902 1 4 Q 1 (4) {2E 1 ) + IOS/Q^Bj)} T 
+ {1052 1 3 Q 1 <3, (2E ] ) -602 1 4 Q 1 (4) (2E 1 ) + 10S 1 5 Q 1 (6) (2E 1 )} <5' 3 
+ { - 1 52/QW (2EJ + 52W (2E,)} § fi 

The coefficients B x , . . . , B 5 can be determined successively by equating the 

'4 ■■ ■ 

coefficients of powers of (5 ; in the expressions for j {Yja)du and — %/5H, If 

«" AC 

additional coefficients B 6 , ... , are desired, they may be found by equating to zero 
the coefficients of powers of ( 5 ; higher than the fifth in the expansion of J (Y/<r) du. 

The expansion of x (R 1? s') in powers of ^ may, of course, be found from the expression 
for Z in the first middle wave without transformation to the Q form. In particular it 
may be proved that B x == 6— 42 1 /cr () . 

24. Second Method of determining the Coefficients J. —When the coefficients B are 
known, the coefficients A 6? . . . may be found in the following way. 

9 t? % 7 

£j JEu m 
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Since x , which is — a (er/er ) 10 3Z/3<r, vanishes at all points of the locus 

r-K, = B 1 ( S -S 1 ) + (B 2 /2 1 )( S -S 1 ) 2 + (B 3 /2 1 3 )( S -S 1 ) 3 + -- ) 
the expression 

945F x (2r) -945o-F/ 1} (2r) +420 

in which 

F 1 (2r) = A + A 1 (2r-2E 1 )+A 2 (2r-2E 1 ) 2 + ... ? 

must become identically zero on substitution of 

Z 1 (B 1 3+B a 3 a +B 3 3 8 +...) 
for r— R x and of 

Z 1 {1+(B 1 + 1)^+B/+B 3 ^+...} 

for <r. Now the powers of cr/Si and (r—'R 1 )lt 1 can all be expanded in powers of $, 
and then the coefficients of the powers of § in the expansion of 

< J4S Lgl - ,«| ™ + «„( -)' ™ _ 105 (|)' ^ 

/V V F 1 (4) (2r) __ (*\* Fi (5) (2r) . 

can be equated severally to zero. The equations thus arising give the values of 
A 6 , A 7 , ... , successively. Suppressing the algebra, which is rather long, we may write 
down the results in the following form : — 

The equation for A 6 is 



2B 1 ^ = 945(B 1 -l)^-210(5B 1 -4)^4 3 + 105(5B 1 -3)2IA 

■Zjj 2-ly 2J\ 2^y 



3 



-30 (56,-2) ^- 4 + 5 (5B!-!) ^T 15 - 

1 * ' 

The equation for A 7 is * 

26^^^ = 26^(945 ^ 3 -945 ^ 3 + 420 4^-105 5 4^ + 15%^- 6 ) 

Z11 \ 2j x Zj X 2j\ 2j x 2j x J 

+ 26^+1) (_945%^ + 840^-315%^ + 60^A'-5%- A 

+ (B 1 + l) 2 (420%^-315?i^2 + 90 4^-10 ^A 

-^ / A 2t A at A 4t A, 5 i A c 6^ A, 
+ BJ 945 |i -1050 4^ 2 + 525 4+^ -150^^+25^^-2^ 

\ c ^— ^1 ^^1 ^^1 ^^1 ^1 ^"l 



6 



4 
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The equation for A 8 is 

4B 1 8 21A = ^(945 ^2-945 %^ 4 +420^^-105^A 6 + 15~ I ) 

■^ „ r ^ %/ 81 A., 4' A, 5' A- fit A. 7'A,\ 

+ 6B 1 2 (B 1 +l) (_945^^+840=^-315 2^+60 2^-5-i-^l) 

\ 2 X Sj 2j 1 2/j 2/i" / 

+ 6B! {B 1 + 1) 2 U20 §^2-315 ^—^ + 90 5iAg-l0 ^-A 6 
+ 3(B 1 + l) 3 (-105^S+60i^-10~ AE 

\ ^1 **i **] 



+ 6B3 2 (945^^-1050^^ + 525^^-150%^ + 25^ 6 -2^4- tN 

\ ^\ ,^1 ***\ "^l '***^1 "^1 / 

+ 3(B 1 +l)B 2 (-1050^4 2 + 1 050^ , -450%^+100^5-10^iA 6, 

6! A 6 \ 

"": < }• 

The equation for A 9 is 



-l A J- -»- 

/ A 9 1 A 3 ! A 1 1 A ' j^ t A 

3B 3 (945|i-1050^^ + 525^^-150^^ + 25~p-2 u - XX6 

\ -Z^i jZ^it JLf-% **4\ *-i\ 



2Bi *9!Aq = 2Bi 4 ^45 41A* _ 94 , 5TA, +42Q 6LA,, _ 1Q5 71A, + lfi 8!A_ 8 V 



8 
/ 



•8 



+ 4B 1 3 (B 1 + l) (-945 ~^ + 840 £1^5-315 %^+60 —A" 5 ^T 

+ 6B/ (B, + 1) 3 (420 %4" 4 - 31 5 ^nr + 90 %r* ~ l ° ^4ir 

\ -wj a^i *-*j "^]_ 

+ 6B 1 (B 1 4-l) 3 -105 ^^ + 60 ^5-10 ^~ le 
4-6 (R 4-1 W 15 — 4 --5 5 

+ l2B/Bj945^4- 3 - 945 ^^ + 420 ^4 5 - 105 %^+ 15 ^T 2 -%r 

\ s x Si Si s x Si 2,1 

+ 6B 1 (3B 1 + 2) bJ-945^ + 840^-315^+60^-5^ 2 
+ 6 (Bj + 1) (3B 1 + 1) B 3 (420 ^2-315 ^T +90 ^5-10 ^ A ° 

+ 3(B 1 + 1) 3 B 2 ('-315^4 S + 180^44-30^~ A 

+ 3B/ (420 ^-525 ^ +300 ^i -100 ^ +20 ^-2 ^A 

/ 91 A q! A -■• 41 A 5' A*- ' "6f\A* 7 ! A 

+ 6B/B 3 945^^^^ 

+ 3(B 1 + l)B 3 (-1050%^ + i0505^-450^ 4 + 100^A-10^ A - 6 ) 



-5 
5' 



sT sj s x 8 • s^ s x * / 

... +3B 4 / 945 A i-1050^ + 525%^- 15 ° i ^ + 25 ^- 2 ^ 

\ Jmf 1 < 4 "^1 "^1 ^^X ^^1 *"fj 

The method avails for the calculation of as many coefficients as may be desired, 
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25. Third Method of determining the Coefficients A. — Another nearly equally effective 
process for finding the coefficients A 6 , A 7 , . . . , is founded upon an expression for t, valid 
at the piston M. 

The equation of motion of the piston shows that at x e = 0, 



du acr 



ii 



dt lOHo-, 10 







and the differential of u is always 



du 7 , du 7 . 
— OOXq-t -777 Cut, 
dXn at 



so that, at x — 0, t can be expressed as a function of s by the equation 

t-T, = - ^1^° f fe)~ n {(Bx-1) +2B^+3B 3 ^ 2 + ... }ds 9 



and thus t— T x can be expanded in powers of S or (s—Si)/^. Also, since t = dZ/du, 
and Tx is the value of £ given by putting r = R x and 5 = S x in the formula for the first 
reflected wave from the left, we have 

,-t, = iQ 5 m-JM)}- Wi m-!a&i 

f F®(2r) F^ (2B,) ! f F^jgr} F" (2RQ 1 f F^ (2r) F^(2RJ 

I o- 7 S 2 7 J 1 <r 6 $!« J 1 o- 5 S t 5 

so that a different form of expansion can be obtained for t—T v By equating coefficients 
of different powers of S in the two forms of expansion we obtain again a series of equa- 
tions giving the values of A 6? A 7? . r . , successively. The results may be recorded as 
follows : — 

The equation for A 6 is 

105|-9(B 1 +l)|i + 2B 1 ^Aj-105J-8(B l +l)^A 2 +2B 1 ^ 

+ 45{-7(B 1+ l)^ + 2B 1 ii^j-lo|-6(B ] + l)^ +2^^} 

+ {-5(B 1+ l)^ + 2B 1 ^} = -10f( / g) 1 °(B 1 -l), 
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The equation for A 7 is 



105 



{-9B 2 +45(B 1 + 1) 2 }|l + 2{ B 2 -9B 1 (B 1 + 1)}^ + | ! B 1 2 ^ 



-105 



+ 45 



10 



{-8B 2 +36(B 1 + l) a }^+2{B 2 -8B 1 (B 1 + l)}^ + | ! B ] 2 iiA.^ 
{-7B 2 +28(B 1 + l)n^+2{B 2 -7B 1 (B 1 +l)}^ + ^B 1 2 ^ 

{-6B 2 +21(B 1 + l)^ilA. 4 + 2{ B 2 -6B 1 (B 1 -il)}^ + ^B 1 26!A( 






V 2! L tf 1 






6! A, 



7! A. 



{-5B 2 +15(B 1 + 1) 2 }^+2{B 2 -5B 1 (B 1 + 1)}^ + |b i ^ 



10 



= -10 



The equation for A 8 is 






o ; 



B 2 -^(B i + l)(B 1 -l)}. 



105 



{-9B 



-105 



+ 45 



i 3 +90(B 1 + l)B 2 -165(B 1 + l) 3 }|i i +2{B 3 -9(2B l + l)B 2 +45B 1 (B 1 +l) 2 }^ 

+ 1 {26^-9^(^ + 1)} ^ + f { B?^* 

{-8B 3 +72(B 1 + l)B 2 -120(B 1 + l) 3 }^A 2 

+ 2{B 3 -8(2B 1 + l)B 2 +36B 1 (B 1 + l) 2 }^ 3 

+ g{2B 1 B 2 -8B 1 2 (B 1 + l)}^ + |B 1 3 ^ 
{-7B a 



53+56(B 1 + l)B 2 -84(B 1 + l) 3 }^ 
+ 2{B 3 -7(2B 1 + 1)B. 



i 2 +28B 1 (B 1 + l) 2 }^ 

2 



-10 



1 
+ |^ {2BA-7B, 2 (B 1 + 1)} ^ + 1^ B/ ^- Al 

{-6B 3 +42(B 1 + 1)B 2 -56(B 1 + 1) 3 }" A 



4! A 4 

?x 6 






+ 2{B 3 -6(2B 1 + 1)B 2 +21B 1 (B 1 + 1) 2 }% A - 5 

1 

+ g {2BA-6B, 2 (B 1 + 1)} ^r- 8 + | B/.^ 
{-5B 3 + 30(B 1 + 1)B 2 -35(B 1 + 1) 3 }^ A - 6 

+ 2{B 3 -5{2B 1 +l)B 2 +15B 1 (B 1 + l) 2 }^ Aj 

+ |{2B 1 B 2 -5B 1 2 (B 1 + l)}^i A -'+|!B 1 3 ^ 

'gj° {B 3 -^(3B 1 + 1) B 2 + 22 (B 1 + l) 2 (B^l)}. 



H/(T Xl ° 

a 
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The equation for A 9 is 



105 



-9B 4 +9.10(B i + l)B 3 +^B/-8 9 - 1 3 ! - 11 (B J + l) 2 B 2 



+ 9 



•10-11 .12 m nV'l^ 



+ 2{b 4 -9(2B 1 + 1)B 3 -9B/ + ^(3B 1 +1)(B 1 + 1)B,- 9 - 1 3 ! - 11 B 1 (B 1 + 1)-J^ 
+ ^|2B 1 B3+B 2 2 -9B 1 (3B 1 + 2)B 2 +^B 1 2 (B 1 + l)j^2 



+ 



+ 1 {8IVB.-9B,' (B 1+ 1)} 114* + J\ Bl ^J 
105[...] + 45[...]-10[...] 

-5B 4 + 5.6(B 1 + l)B 3 +^B/-3^| r ^(B, + l) 2 B 2 + 5 - 6 4! 7 - 8 (B i + l)j^ 

+ 2{b 4 -5(2B 1 + 1)B 3 -5B/+^(3Bh-1)(B 1 + 1)B 2 -^^B 1 (B 1 + 1)4^4 ? 
+ |^{2B 1 B 3 +B/-5B 1 (3B 1 + 2)B 2 +^B 1 2 (B 1 +1) 2 }^ 

+ f { {SB/B.-SB, 3 (B l +l)}?14 s + ^B 1 * ^ 



= -io!fe xl ° 

a\Z v 



B 4 -^(2B 1 + l)B 3 -^B/+i^B 1 (B 1 + l)B : 



11.12. 13 (Bi + 1)*/-d .v' 

— si — (Bi-1) . 

where the law of formation of the terms that are not written down is sufficiently obvious. 

The formulae of this article may, of course, be transformed into those of the previous 
article by means of the relations by which the coefficients B were expressed in terms 
of the differential coefficients of Q 1? and the relations by which the coefficients A l9 ... y A 3 
were expressed in terms of the same differential coefficients. They are useful in 
numerical work as affording a verification of the values obtained for the coefficients 
A. ? A 7> . . . , from the previous formulae. 

Formulae similar to those of the present and preceding articles may be obtained for 
the coefficients in the expression for Z belonging to the first reflected wave from the 
right, but it seems hardly worth while to write them down. 

The Second Middle Wave. 

26. Method of determining the Second Middle Wave.- — The first reflected wave from 
the left meets that from the right at the place and time determined by substituting 
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E x for r and s t for s in the formulae belonging to the first middle wave. When they 
meet, the first middle wave becomes obliterated, and the second middle wave begins 
to be generated at the time and place in question, and encroaches upon the two first 
reflected waves. 

To determine the second middle wave we have the conditions that at its advancing 
front, where r = R 1? the Z belonging to it is equal to that belonging to the first reflected 
wave from the right, and at its receding front, where s = s u the Z belonging to it is 
equal to that belonging to the first reflected wave from the left. Riemann's method 
may be applied in exactly the same way as in Article 15. If P is the point (r\ $'), A 
the point (R 1? s'), B the point (R l5 s^ and C the point (r 7 , s 2 ), we have 



Z(r\ s') = [YZ] A -[YZ] B + [VZ] C + f Z 

Jab 



'3V 



ds r + sj 



5Y )ds-\ a 



'3V 



5V\ 



BC 



dr r + sj 



dr. 



At A we have 



r = H 



i? 



s = s , 



f = o, 



TT ___ / -tw] I & 



r'-fV 



and Z is the result of substituting R 2 for r and s' for s in the formula 



Z 



1 3 X4 

/^ + l x U + ( — 



cr ocr, 



At B we have 



CCT, 



10 



3o7 Lcrl945 



a 



+/i (° r -^) 



r = R 1? * = *, g = - ffi- nUV't V = (^±^Y(l-20f+90?- 140^ + 70^), 

(r + 5') (li! + 5J \ r+s J 

and Z is the result of substituting R x for r and s x for s in the formulae for the first middle 
wave, or in those for either of the first reflected waves. For the present we shall denote 
it by Z B , and observe that it is independent of r f and s'. 
At C we have 

r = r\ s = s u f = 0, Y = 



V + s^ 5 



y + $' 



and Z is the result of substituting r f for r and s x for s in the formula 



Z = Jii + J-ii^4" 



1 9 V fFjjfor + M)') 



cr 3cr, 



O" 



Along AB, where r — Rj and s increases from s' to s u we have 



Z = yfc 1 +z 1 (R 1 — «)+ 



i a 



,CT '3<T, 



1 J£2jl1 +/ (o._ w ) 
.cr 1.945a ,/lv 



3J_ 5V = (R 1 + /)(R l -/)(R 1 +^ (20 _ 180 ^ + 420fj280f)) 



os r + s 



s = 



(r' + s'Y 

(B,-r')(s-s') 
(r'+s')(R,+s)' 



VOL. OOXXII. — A. 



2 W 



200 MESSRS. A. E. H. LOVE AND P. B. PIDDUCK 

and along BC, where s — s x and r decreases from Rj tor', we have 

|^-— = (^ + r/ )(^-^)(gi+^) ,i (20-180^+420f-280f), 
or r + s (r + s) 

» (r-rQfa-sQ . 

(r' + s')(r+s 1 ) 

The value of Z at (r', s') can be regarded as a sum of terms with the coefficients 

b> *^i> '1? ^^"o fi'Q*y(X / ) ^oj ^i) ••• ; -*^-i) -^l) -^0) -^i) ••• » 

and each of these terms may be found from the formula for Z (r\ s f ) by performing the 
integrations where necessary. The result will be to exhibit Z (r\ s') as a sum of terms 
with these coefficients. 

27. Determination of the Second Middle Wave. — No integration is needed in order 
to obtain the term which has Z B as a factor, but it is important to observe that V B , 
as a function of r' and s\ can be expressed either in the form 



v B=l(Bl+8l) (^y{ (*--s.>y + «.> 



4 



or in the form 

v,=KB. + *,AArf (/ - B - ) y +s -n - 



a 



<T OCT/ I cr 



We shall suppress the accents on r ' and s' so as to express the value of Z at (r, s). The 
term with coefficient Z B is 



The term with coefficient h x is 



p/IAYj(£z^H£±M 5 j. 



cr do- 



er 

The term with coefficient l x is 



' 51 W 3cr/ 1 cr J 



The terms with coefficients ccr 10 /945a and a are 



. ^ + «.) a 3 



945a / \<r 3cr, 



1 / s— s* y L , , / s + 



•1 



1+4 



a VKi + v I ' \E 



,! + V \E X + 8 J 



VRj+V \Ri+*i/J_ 
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The terms with coefficients a u a 2 , a 3 are 



— 2a x (R1 + S1) 



i a 



.cr 3o\ 



i / s-s 1 y\ 1 1 g ^ + R^ , 6 /g+Bi 



_cr VRj + S] 



Ri + ^Z ' "VRj + ^i 

R\ o 
. 1 



+ 10 



2 2 a 2 (E 1 + s 1 ) i 



i a 



2 



1( *-»iYl 1 + 2 (°^) + 9 (°^)' + i(p*L) +5 



cr3o-y Lo-VEx + g, 



Ei+«i/ VEj+Sj, 



■ 2 8 a,(R 1 + * 1 ) 3 (i^ N * 

\cr do 



1 / S — S 1 



.cr\R l + S lt 



1-f 



s+RA , /s + RA 2 , /s+Rj 



+ 



+ 



S + Rj 

vRi+^jL 

s + R^ 



.Ri + 5i 



+ 



3i+5i/ \Ri+v xRi+v VRi+^i 



r 



The terms with coefficients a 4 , a 5? . . . , are 



4/I 8V[ (^^ 

4 \cr da \ (T 



2*a 



2°a, - 



1 8 y [ (a'-fr) 1 

s cr 3cr/ \ cr 



and so on. 

The term with coefficient Ki is 



fKi( iay { (^^n. 



The term with coefficient L 2 is 



_-_ ALi ( ± Jl ) i (r-R 1 ) 4 (5^ 1 -4R 3 ^r)(r+6- 1 ) 5 ] 



1 _a X4 

cr 3o\ 



cr 



The terms with coefficients A , A u A 2> A 3 are 



A 







1 avri/r-EAU 1+4 ^z±£LUio('^±^Y + 2o( / ^±^y + 35^ +s ' 



cr 807 [_a\Ri + $: 



— 2Ai(Ri+^i) 



1 a 



XT OCT, 



"1/ r-Bi V 

.cr \Rj + 5i, 



1 + 3 






+ 6 



VRi + 5 1( 



Rl + Sv 



JRi+si/ \Ri+5i 






2 2 A 2 (R 1 + 5l ) 2 (^^ 4 

\cr 8<x 



1 /r-RiV 



2»A 3 (R 1 + s 1 )< 



cr \R L + Si 

1 a X4 



1 + 21 



r + Si 



+ 3 



/ r + Sl \a 



,Ri + si/ \Ri + «i, 



+ 4 



>Ri + *i, 



+ 5 : 



r + ^i 



41 



Ri + sj J _ 



cr 3cr y 



l/r-E 1 Vi 1+ /r ±£L 



r + 6'i X2 



.cr VRj + 5i, 






^ + ^i V . / r + Si X4 



H + 



R1 + S1/ \Ri + 5i/ \Ri + 5i/ • \Ri + Si/J_ 



The terms with coefficients A 4? A 5 , . . . , are 



\<r cor I cr 



4 



6 \o-8cr/ 1 o- 



and so "on. 
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28. Indication of a General Method. — If the coefficients a or A with suffixes exceeding 
10 were all zero, the expression for Z in the second middle wave could be transformed 

at once to the form 

(I 3 V fQ 2 (er+^| +K +L ^ 
\cr ccr/ I cr J 

or to the form 

G" OCT J t <T J 

But if these coefficients do not vanish a transformation of the same kind is still possible . 
We have, for example, 

1 gyf (o-+^)" + (a-< | = 2 5 5! ((j2+lM 

,cr3a7 I <r J 

(ld\U<r + ur-{*-uT\ = 275 , u(3(r * +nu% 

\or dor J [ cr j 

and thus Z can be expressed either in the form 

^j(S^)} + K, + L I „ +H ,(»= + n»-) + N,»( 3 ^n»') + ..., 

cr ocr/ I cr J 

or in the form 

_ _\ < 23l_ri_il +^ 2 4.^4. m3 (o- 3 +ii^ 2 ) 4. n 2 u(3or 2 +llu 2 ) + ..., 

where the factors whose coefficients are written K 2? L 2 , . . . , or Jc 2 , l 2 , . . . , are the 
homogeneous rational integral functions of cr and t£ of degrees 0, 1, 2, 3, ... , which 
satisfy the differential equation for Z. 

When this transformation is effected we may proceed to determine the second reflected 
waves. The first step is to find sets of coefficients analogous to A , ... , A 5 and 
a , ... , a b . The next is to find sets of coefficients analogous to B 1? B 2 , ... , determining 
the loci in the plane of (r , s) along which x = and x = c during the time that these 
reflected waves are being generated. By means of the coefficients analogous to B l5 
B 2 , ... , sets of coefficients analogous to A 6 , A 7 , ... , and a 69 a 79 ... , may be found, and 
thus the second reflected waves may be determined. 

From the formulae for Z in the second reflected waves that in the third middle wave 
may be found, in the same way as the formula for Z in the second middle wave was 
found from those in the first reflected waves. 

The method of solution can be continued, and gives a theoretically complete solution 
of the problem ; but when arithmetical computation is attempted, failure may arise 
through approximate equality of groups of terms with opposite signs, so that some 
quantity, which ought to be calculable to five figures, for example, may only be calcu- 
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lable to three. This difficulty was found to present itself in the calculation of the second 
reflected waves by this method, and another method had to be sought. An account 
of this will be given in the theory of the second reflected waves. 

29. Pistons of Equal Mass. — A considerable reduction in the number of coefficients 
to be calculated is effected by supposing the two pistons to have the same mass. When 
this is so H = h, and hereafter we shall write everywhere h for H. The calculation of 
A , A l5 ... , A 5 is then simplified a good deal. Further, it appears that the coefficients 
a differ only in sign from the coefficients A, or we have 

It is now unnecessary to calculate separately the pressures, velocities, displacements, 
and times at the two pistons. We shall speak of the piston specified by x = c as the 
" shot," and of the piston specified by x = as the " image of the shot." We shall 
generally calculate the pressures, &c, for the image of the shot, because a slight simpli- 
fication is effected by putting x equal to zero. 

30. Incidence of the Second Middle Wave upon the Pistons. — The value of s at the 
receding front of the second middle wave is that which has been denoted by s u and 
in the case of equal pistons it is the same as Rj or Si — Jcr . This is therefore the value 
of s at the image of the shot at the instant when the receding front of the second middle 
wave reaches it. It will be denoted by S 2 . The corresponding value of r may be found 
from the formula 

r-B, = B 1 («-S 1 ) + (BA)(«-S 1 )» + ... 

by putting S 2 for s. It will be denoted by R 2 . From this the corresponding value of 
or may be found. It will be denoted by 2 2 . The corresponding value of u 9 which is 
R 2 — S 2 , will be denoted by U 2 . The corresponding value of Z, denoted by Z 2 , can be 
found most simply from the formula for the first reflected wave from the left. We 
have 

Z, - K 1 + L 1 U, + 1 «»5 ^ -105 mi + 45 EqfflJ-io ™> t B|M, 

^2 ^2 ^2 ^2 ^2 

where 

Fi (2B 2 ) ' /%\ w J A, _ Ai / 2Eg - 2B 2 
Z™ W IS, 10 2A 2i 

1 2 ! A 2 / 2Bi - 2R 2 Y 1 3 ! A, / 2B t - 2B 2 V , , 
+ 2! Zf { 2, / 3! ^ V Sr J + -J' 

F 1 ( »(2R 2 ) = /2A 9 J Aj _ 2JA, / 2Ri-2B 2 \ 

£ 2 « WW ^ \ ^ / 

1 3 ! A 3 / 2Bx - 2B 2 \ 2 1 4 ! A 4 / 2Bi - 2B_ 2 \ 3 , 1 

+ 2\~$r\ % ) 3! v { $t / + -j ' 
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F 1 (a) (2R 2 )_/Sr\ 8 /2!A ! 



S, 8 



V 



3! A 3 /2Hi-2E 



F 1 (3) (2E ii ) _ /SA 7 f 3 ! A : 



Zi 7 \ ti 

1 4!A 4 / 2Ex-2 E 2 \ 2 
+ 2! ^ \ Si 

4! A 4 /2R!-2E; 



1 5!A 5 /2E!-2E. 



3! V 



Si 



-f" . . . r , 



2, 7 



^5»/ L ^1 



2> 6 



^1 



+ 



1 5!A 5 /2Ei-2E 2 \ 2 1 6IAJ2EX-2EA 3 



F! (4) (2E 2 ) 



/SA 9 |4!_A 4 

\^2' v. ^1 



2! V \ S x 
5! A 8 /2E 1 -2E 8 ' 



3! ^ 



Si 



+ — r, 



^ 



Sx 



+ 



I 6!A e / 2E 1 -2E 2 \ 2 _ 1 7 ! A 7 / / 2E 1 -2E 2 \ 3 



2! Sx ! 



Si / 3! V 



2*i 



+ ...K 



31. Transformation of the Formula for the Second Middle Wave. — In what follows 
we shall disregard coefficients A beyond A 9 ; if it were desired to include further co- 
efficients A some of the formulae would require modification but there is no difficulty 
arising from the convention to stop at A 9 . The most effective transformation of the 
formula for Z in the second middle wave is found by putting for Z B the value derived 
from the first reflected wave from the right, viz. : — 



Z B = £i + ^i(Ri-si) + 



1 3 V [ co- l ° a + a 1 (2s-2s 1 )+a 2 (2s-2s 1 ) 2 +... 



o" 3 



<rj l945acr 



+ 



<T 



r = Ei, 8 = Si 



so that the terms contributed to Z by Z B come to 



f UiORi+^ + ios 



co-, 



10 



945a 



+ a ) —105 (R x + Si)ai 



x 



+ 45 (R^*) 2 2 ! a 2 - 10 (Ri + Si) 3 3 ! a 3 + (Ri + Si) 4 4 ! a, 

i a v r (g-gi) 4 (g+Ri) 4 l 

crCRi + ^y 



cr3 



cr, 



k8 



and then, before putting Ei for s u or ■— A , — Ai, ..., for a 05 a 1? ..., transforming the 
terms contributed by A , A ]? A 25 A 3 to the form 



cr do-. 



i {A0 + 2A! (r-R x ) + 2 2 A 2 (r-R 1 ) 2 +2 3 A 3 (r-Rx) 3 } 



+ - 



IJLY 

cr 3<t/ 



cr 

"A 



56 



Si + r 



Ri + Si, 
2Ai(Ri + Si) 



+ 140 



s*+r\ 6 10A / 6'! + r V + 35 / «i+r 



'8' 



\JLi)i + <Si 

5 



-120 



CT 






2 2 A 2 (R 1 + ^ 



■21 ( 



6 



/ Si+^ \° , xe I Si + r 



Ri + Sl 

6 



Ri + ^b 

7 



Ri + «l 



+ 56 



Jtti + 6*i 



50 



Sx + r V | j5 ( 6 'i + r 

\jL\/x + #1, 



8" 



cr 



*i + r\ 5 +17 ^ : ji + ry -16 



JEti + Si, 



2 3 A 3 (R 1 + s 1 ) 3 [ / ^ + r . +3 

VRi + V VR X + i 



5 / Sj + r 



Ri + s v 

6 






3 



5i + ? 



,Ri + s lt 

Si + r 
Ri + Si 

7 / Sl + r 



+ 5 



+ 



S irti + Si/ \JrCi + 5i 



5i+r 
,Ri + Si 

8" 

I. 
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The first line of this expression, with the terms contributed by A 4 , A s , . . . , makes up 

\cr da) 1 a 

and the remaining lines are unaltered when — s is written for r. 
The terms contributed by K* and Lj are the same as 

Kl+LlM _ l K,(iAy{(^zia±ii] +ALl (iAY{(£z£ ! m ! i!i 

\a da) I or J \a dor J I a 



5" 



and thus the Z of the second middle wave is expressed entirely as the sum of the Z of 
the first reflected wave from the left and a function of the form 

K a da) \ or 

Further, noting that with equal pistons ^ = L l5 we see, that <p 2 (a--u) contains no 
terms of degree higher than the ninth in s or ^{a—u)l Also we see that it can be 
expressed as a rational integral function of (s— Si)/(Ri+Si) of the ninth degree, and that 
it contains no terms of degree lower than the fourth. Since Z and dZjdo- are continuous 
at s = Si with the Z and dZfda belonging to the first reflected wave from the left, the 
function <p 2 can contain no terms of the fourth or fifth degree in (s— Si)/(Ri+$i). The 
vanishing of the coefficients of these terms does not introduce any new condition. On 
replacing a , ... by — A 0? ... , we have the result that in the second middle wave 

where £ is written for (s— Si)/(Ri + Si), and s 1 = Ri = .'S a , while j/ 6 , w, %, % are given by 
the equations 

\ 6 = -|(K 1 -^^^-^(| ? Y U + 280f + 140^ + 62^ + 23r3+6^-f 6 , 
w = 4(K I -* 1 ) ^!-^(gJ°-240£ ) -120f 1 -52&-18&-4&-&, 
% = -KKx-fe) (^L _-£_/|>Y + 70£ ) + 35f 1 + 15&+5fi, + ^-£ J> 



w = |(K 1 -* 1 )^. 9 -& > 
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in which 



? _ A fa) 10 t _ 9 A X 2K, faY t _ 9 2 ! A 2 /2RA 2 fa 



t — — 4 

S3 — — i 



4 Sj L A 3 (2R } Yfa\ 1Q , _ 2 jyu /gBAVgiV * _ _ 4 5 ! A 5 /2RA VSA 10 

* sj \ ^ / W ' i4 " 3 Si 6 \ ^ / W ' ^ 5 " T * s/ \ t x 1 \%J 

, „ 4 6!A 6 /2RA 6 /SA 10 , _ 8 _7!A 7 /2E 1 \YS 1 Y° , _ __2 8! As^RAY^A 10 

U - 45-^r-^y ^J , fir - -3T5 -tf\-tj \ZJ > & - ^^^T\X/ \%) ' 

t _ _4_ 9 ! A 9 /2Ex Y /SA 10 



The Second Reflected Waves. 

32. Relation between Pressure and Velocity at a Piston. — The relation between pressure 
and velocity at the image of the shot is an equation connecting r and s, which holds at 
x = 0, and can be interpreted as the equation of a certain locus in the plane of r and s. 
This equation can be written in the form 

(r-R 3 )/2 2 - B' 1 J + B'^ + BV+ 

where § stands for (s— S 2 )/2 2 , and the coefficients B' are at present undetermined. 

To determine these coefficients we have recourse to the method of Articles 21-23. 
During the progress of the second reflected wave from the left, the value of x at any 
point in the region occupied by it can be expressed in terms of the values r r and s r of 
r and s, which occur simultaneously at the point, by the formula 



x = — 5 A 



AC 



(Y/cr) du, 



wherein the integral is taken along the locus from the point A, where r = r' 9 to the 
point C, where s = s'. In this integral 

u-JJ, = $ 2 {(B' l -l)S+B' 2 ?+BW+ ...}, 
du = t 2 {(B'i-1) +2BV+3B / 3 5 a + ...}<M, 
o- = S s {l + (B / 1 + l)5+B' s #+By 8 + ...}, 

I^i^^J 1+8 o(70(^ +210 (7^(;-f)' + ... 
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On putting R 2 for r', we have the value of x along the junction of the second middle 
wave and the second reflected wave from the left expanded in a series of powers of #, 
or (s'— S 2 )/2 2 , in the form 

xo = -5h[iB\-l) $' + {B',-(3B'i-2) (B'i-1)} #> 

+ {B , 3 -3B / 2 (2B / 1 -l) + (7B?-6B / 1 +2)(B / 1 -l)}r 
+ {B' 4 -2B ; 3 (3B'!- l)-3B' 2 2 -|B' 2 (42BV-33B' 1 + 5) 

-J(28B?-2lBV+9B / 1 -2) (B'i-1)} S fi 
+ {B' e - JB / 4 (30B' 1 -7)-6B / 2 B , 3 + t 1 oB / 3 (210B?-108B' 1 + 7) 

+ ^B7(l05B' 1 -27)- T VB / 2 (560B / 1 3 -462B?+8lB' 1 -7) 

+i(l26BV-56B' 1 3 +2lB7-6B / 1 + l)(B / 1 -l)}^ 5 + ...]. 

Now at any point (R 2 , s') on the same junction the value of x can be obtained by 
forming — II 3Z/3cr, where 

and 

n = a (cr/o-o) 10 , i =(5-S a )/2E! f 

and putting therein 

cr = B 2 +s', w = R 2 — s', 5 = 5', s'—S* = 2 2 <5', 
and the result can be expressed in terms of <$' in the form 



Xo _ (%\ 



^ . Vo 



10 



. I 2 2 . ^2 ^2 

_ 105 (l + a7?^i^+^^ 

Z 2 ^2 ^ 2 

+ %(^-) 6 {945^-945 x 3 (l +§') <T 5 + 420 x- 1 / (l + 3 , )"3' 4 

- 105 x 15 (1 + <S') 3 <T+ 15 x -V- (l +.S') 1 ^-^ (l +<?? <*'} 

- w ( / 7 %-Y{945^ 7 -945x|(l + ^)r + 420x-^(l + <5 / )^ /6 

_105xJ^(l + ^) 8 r+15xJ^(H-^) 4 r-H i (l + ^) 6 ^'} 

+ /_|_Y {945^-945 x 4 (l +<$') <T + 420 x 14 (l + ^) 2 ^ 

-105x42(l+^) 3 ^ /5 +15xl05(l + ^) 4 ^-210 (l + ^') 5 H 



m 



A 2 \ r^.rt/l ^> i r « / 1 . */\ W8 . i(m..lo/l i }'\2 « 



{945a /9 -945 xf (1 + ^) r+420x 18 (1 + ^' 



2Ri, 

-105x63(l+n^ /6 +15xl89(l + ^)^ /3 - A f-(l + ^) 5 «5 /4 } 
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The terms independent of $' in the right-hand member of this equation add up to zero, 
for x vanishes at (R 2 , S 2 ) ; and, by equating the coefficients of powers of $' in the two 
expressions for x 09 equations are obtained from which the values of the coefficients B / 
can be found successively. 

33. Relation between Velocity and Time at a Piston. — The time at which any par- 
ticular simultaneous values of r and s occur at x = 0, can be found by the method of 
Article 25, and thus the relation between velocity and time at the image of the shot 
may be traced. We can write down the equation 

t-T 2 ■= -10* |Vi f fey il {(B / 1 -l) + 2BV+3B / /+ ... }ds 9 

CO 2<2 J Si ^3/ 

in which 

S = {s-8*)p,2, c = 2 2 {l+(B / 1 + l)^+BV 2 +By 3 +... 

and thus £— T 2 can be expanded in powers of S in the form 



where 



lOh 



CTc 



10 



•^-—(rYV.-v-w+iMB'.-i)}, 

(JO \^2 / 

o 3 = - ^fe Y° {B' 8 - V- (3B' 1+ 1) B' 3 +22 (B' 1 + l) a (B'i-1)}, 



a \2 



3 



o 4 = - ^fe ) 10 {B' 4 - V- (2B' 1+ 1) B'.-JgW+eeB'x (B' 1+ 1) B' s 



_xfii(B' 1 + l) 3 (B' 1 -l)} ) 



10A 



(To 



10 



c 5 = -^(1° ) {B , 6 -^(5B' 1 + 3)B / 4 -llB' 2 B' 3 +J T 3 T 2 ( 5 B' 1 + 1) (B', + 1) B\ 



3 
2 



+ W (5B' 1 + 3) BV--W (SB'a- 1) (B' 1+ 1) 2 B' s 



34. Displacement of a Piston. — To obtain the displacement of the image of the shot, 
we have to find the value of x at x = in terms of simultaneous values of r and s 
occurring on the locus 

r-E 2 = B\ (*-S,) + (BVS,) {s-8,)*+{B'Jt 2 2 ) (*-S a ) 8 + ... . 

Now, when # = 0, we have x = ut—Z, and for the value of x at (R 2 , S 2 ), denoted by 
X 2 , we have X 2 = U 3 T 2 — Z 2 , so that when x = 0, we have 

^-X 2 = ^^U 2 T 2 ~(Z~Z 2 ). 
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Also we have 

Z-Z 2 = 



%} 



^ az , , az , 

(r 2 ,s 3 )w do- 



where the integral is taken along the locus, and, since 3Z/3<r vanishes along the locus, 
we get 

Hr t s) rs flf 

__. npi i 7 1 1*1' 7 

J(R 2f S 2 ) Js 2 rf^. 

in which 

^ = U 2 +S 2 {(B , 1 -1)(5+B / / + BV 3 + ... }, 

$ = T 2 + Ci(5+<vP + c 3 <5 3 + ... , 

It follows that, at # = 0, x is given by the equation 

x-X 2 = IW+ {U 2 c 3 +^ 2 c 1 (B^l)} #+ {U A +^ s (B , 1 -l)+iS a cB , jr } <5 3 . 
+ {U 3 e 4 +t2 2 c 3 (B / 1 -l) +i^B , 3 +it2C 1 B / 3 }(5 4 
+ {U2C5 + £ S 2 c 4 (B'i — 1 ) + f S 2 ^sB / 2 + f ^^B's + iS 2 ciB / 4 } (f + . . . , 

which may be written 

The formulae which have been obtained avail to determine the displacement, velocity 
and pressure at the shot or its image at any time during the progress of the second 
reflected waves. 

35. General Method for the Second Reflected Waves. — We shall need to be able to calcu- 
late Z, t and x for any simultaneous values of r and s that can occur in the second 
reflected wave from the left. It is best to obtain formulae for t and x separately, and 
not to deduce them from the formula for Z by differentiation, because the formulae 
will be approximate, and to obtain the terms of any particular order in t, for example, 
by differentiation of Z it would be necessary to obtain the terms of order higher by one 
in the formula for Z. The method of determining the formula for x has been indicated 
already in Article 32, and the work will be completed presently. The formulae for t 
and Z will be found by similar applications of the method of Eiemann. We begin with 
the formula for t. After finding formulae for x 09 1, Z, we can calculate x from the equation 

—*{! + (?) , ( 1 -!)H- z - 

36. Method of Determining t.— The value of t along the locus x = has been found 
in Article 33. To obtain the differential coefficients of t along the same locus we have 
the equations 

uXq yj ut __ ^ vXq yj vt s\ 

U dor dcr du 

Ji Gr Zi 
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and 

in/ ^1 

dx „ 10h ? Sx du 

dcr f -. /<^or*\ 2 ] du 



m\ m *nm 



where do-jdu is to be found from the equation of the locus. Thus we may write 



d 



cr 



dt = her 1 " du 3* _ 1ft W° 1 

3<r coo- 11 ~J^\ 2 ' du acr 11 1 _(d^) 2 ' 

\du ) \du I 

Now let (r\ s') be any point P, and let lines PA, PC parallel to the axes of s and r 
meet the locus in A, C, as in fig. 3 in Article 22. Then, since t satisfies the same differen- 
tial equation as Z, the integral 

j \cs cr/ \cr cr/ 

taken round the contour formed by the arc AC and the lines CP, PA vanishes, and 
therefore we have the equation 

t(r\ s>) = [Yt] k+ f (|^-5 ¥) t dr + (|< +5 i) Vcfe, 

J AC \ or <TJ \OS (T/ 

or, on putting 

the equation 

*' (r>, s>) = [ W] + f ( 8 J--5^) t> dr + (|*-' + 5 t -)v ds. 

But we have, along the locus, 

dt' _ dt' dt' W° 1 _ 5W° du 



' "i IMtlMUMMWiBlM 



3s da du acr 11 dcr 1 acr 11 ds 

du 

and, by the theory of Article 25, the expression last written is the same as the value of 
^dtjds along the locus, or we have 



Also along the locus we have 



* / = c 1 a+cJ 8 +c«3 8 + 
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Let 5 A denote the value of s at A. Then at A we have 



'A 



r = r, s=:s A , Y = {^r^j, 

where S A stands for (s A — S 2 )/2 2 . 
Along AC we have the formulae already written for t' and dt'jds, and we have further 

\r + s J (r + s) (r + s) 

§V 5 Y = ^ + r/ H^^i r + ^^ 20^180^+420r-280^ 

and we have to put 

<fc = $ 2 d8, dr = S a (B / 1 + 2BV+3B' 3 3 8 +.;.)^ 

The limits of integration are 5 A and ^, which is ($'— S 2 )/2 2 « 
The value of <5 A is to be found by reversing the series 

r-R 2 = (s-S 2 ) {B' 1+ (B' 2 /2 2 ) (s-S 2 ) + (B' 8 /2/) (s-S 2 ) 2 + ...}, 
and putting r f for r and s A for s. If we write / for (r'—R 2 )/2 a the result is 



& _ _£ JJ 2 /2 , I jjj 2 __ ^> 3 I '3_/ 5-P 2 __ 5j3 2 1> 3 , L> 4 

A " B\ BV e 1 BV BVJ e \ BV B? BV/ e 



/ TV 4 TV 2 R' "R' 2 TV TV TV \ 

j. Mi ±1JL_ 21 2 3 1 a 5 2 I fi 2 4 5 1 j/ 5 I 

V B ; 9 B' 8 B' 7 B' 7 B' «y 

Thus S A is known in terms of r f . 

37, Formula for the Time. — We work out the formula for t' or t— T a in terms of $, 
or (s — S 2 )/S 2 , and e, or (r— R 2 )/£ 2 , at any point answering to simultaneous values of 
r and s which can occur in the second reflected wave from the left. For this we first 
perform the integrations with respect to r and s and then suppress the accents on cY 
and e. We record the results as far as terms of the fourth order* 

The terms of the first order present themselves in the form 

and it is simpler to leave the factors 

/ l+e + 4 V 
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as they are, rather than to expand them in powers of e, S A and S. In like manner the 
terms of the second order are 

lT^tf cA2+ (u^W { ^ + ^ (B ' 1+8)} ( ^ 3) ' 

the terms of the third order are 

T^if eAS+ (i+.+j)' { * Cs+ |c * (Fi+2) +hi (B ' i+ 1} (B ' i+ 5) + |CiB ' 2} (*-*A 

+ ^rw Gie {S ~ SJ " (rrbr ClFl {s+ 2 ' a) ( ^ a)2> 

and the terms of the fourth order are 

TT^ff cA * + (i+l+sy ^ + *^ ( 3B/ *+ 5 ) +fa ( B '*+ *) ( B/ i+ 3 ) +fo B '* 

+f c x (B'x + 1) 2 (B'x + 7) +|* (B' a + 3) B'^+foB',} (<S 4 - V) 
1 {W^+|*B\(5B\4ai)+^^ 



(l+e + 5) 



38. Formula for Z.— -We write Z' for.Z— Z 2 , and seek first a formula for Z' along the 
locus x == 0. The value of Z' along this locus is given by the equation 



C(r,s) g2 3Z 

Z' = — -dcr+ — du, 



where the integral is taken along the locus 3Z/3cr = 0, so that 

Z'.= f t^ds = sJ S (T 2 + c 1 (5+c/+...)(B / 1 -l + 2BV + 3BV 2 -f ...)(«. 

Thus the value of Z' along the locus can be expanded in powers of S in the form 

Z' = ^5+^+^+..., 



where 



4 = $ 2 {T^+fciB'.+ic^B'i-l)}, 

d 4 = 2 2 {T^+fcxB's+lc^+i^ (B'i-1)}, 

d 6 = ^{T^+lcjB^+lc^+lcaB^+lc^B^-l)}, 
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We require also the differential coefficient dZ'/ds along the locus, and this is given 
by the equation 

— - = — t = — (l 2 + Cio + C 2 o + ...). 

The value of 71 at any point (r' 9 s f ) is then given by the equation 

71 {r\ s f ) = [VZ ; ] A + f f f^ - ~ ) Z' dr + (^ + ^) V tfo, 

J ac \ or cr / \os or j 

in which V has the same form as in Article 36, the value of 71 at A is 

and the integration is taken along the locus. 

The result may be recorded in a similar form to that for t f in Article 37. The terms 
of the first order in the formula for Z' are 

___|j <y A _ (l + e+5)8 S a T a (*- J A ), 
the terms of the second order are 

the terms of the third order are 



TX-Z7 ^ 3 - rr - 1 nv 5 [iS a {c. + 5c 1 (B^ + l) + 10T 2 (B^ + l) 2 +5T 2 B^} 

1+e + d/ 11+eH-o) 

-lH+4d 1 (B' 1 + l)}](^-4 8 ) 
+ (i +!- f .^ ¥ ' (^ T »-^> B '* (*+ 2 **) ^~^) 2 " (i +«+ j)' ^ ( ^~^ a)2 ' 

and the terms of the fourth order are 

+ 10T 3 (B'x + 1) 3 + 20T a (B'i + 1) B' a + 5T 2 B' 3 } 

- i{d s +4d a (B'x + 1) + 6^ (B'x+ l) 2 +4d 3 B' 2 }] (J 4 -4 4 ) 

[¥^{oi + 4T a (B / x + l)}+-W 1 (B'i-5)]e(^+24)(^-0 2 



I J[ _|_ £ _j- {} I 



+ /, 1 n« [P 2 {cxB' 1 + 4T a B'x (B'x + 1) + T,B' a } -f {d,B' 1 +*B' 1 -(3B' 1 + 9) + 2^,}] 

x($*+2§J+3S A i )(8-S x y. 
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39. Formula for x . — A formula has been obtained in Article 32, and can be written 
in the form 

Xo ( r ',s') = - 5h f /rl±^Y j i + 3 o (^^1%=!^) + ...}^(B / 1 -1 + 2BV-F-3BV 2 + ...)dS, 

J5 A \ o- / L \r + s) a- J a 

in which we have to put 

■ r-V = Z 2 {B\S-e' + B' 2 S 2 + B' 3 S 3 + ...}, 

After the integrations are performed the accents on $' and e are to be suppressed. The 
result may be recorded in the form : — -The terms of the first order in the expression for 

-5h{l+ e +Sf{B\-l){S-S A ), 

the terms of the second order are 

-5h(l + e+SY{B' 2 -3(B\ + l){B' 1 -l)}{S 2 -§ K 2 ), 

the terms of the third order are 

-5^(l+e+^) 5 {B / 3-2(3B / 1 + l)B / 2 +7(B' 1 + l) 2 (B / 1 -])}(^-^ A 3 ) 
-150h (1 + e +SY {£ (B' x -1) e-}B\ (B',-1) (8 + 2S*)} (S-SJ, 

and the terms of the fourth order are 

-5h(l + e+§f{B' i -S{2B\ + l)B , s -SB' 2 2 + 2lB\{B' l +l)B' 2 

-14(B' 1 + i) 8 (B' 1 -l)}(5 1 -V) 

-25A(l+e + ^) 4 [{2B' 2 -7(B' 1 + l)(B' 1 -l)} e ( ( !+2^) 

-|{(3B' 1 -l)B / 2 -7B' 1 (B' 1 + l)(B' 1 -l)}(«5 2 + 2^ + 34 2 )](^-<5 A ) 3 . 

40. $£ate of <Ae (?a^ a£ any Time.— With a view to applications it is important to 
indicate how the state of the gas may be determined at any time, or when the shot and 
its image have both travelled an assigned distance. We shall suppose that the time 
in question is an instant during the generation of the second reflected waves, before the 
second middle wave is obliterated. Then the central part of the tube is occupied by 
the second middle wave, and beyond the junctions the rest of the tube, up to the shot 
and its image, are occupied by the second reflected waves. 

An assigned position of the shot and its image answer to a given value of x, and the 
corresponding value of S is to be found by solving the equation for x given in Article 34. 
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This is the value of S at the image of the shot, and the corresponding values of r and s 
at the image of the shot are given by formulae in the same article. Also, $ being known 
for the image of the shot in this position, the value of t is given by the formula of Article 
33. Let this particular value of t be denoted by T 3 , and in like manner let the values 
of the various quantities at the image of the shot at this time be denoted by attaching 
a suffix 3 to the letters, thus : — R 3 , $ s . 

In the second reflected wave from the left the values of r that occur lie between E 2 
and B 3 . To each such value, when t = T 3 , there answers a value of s and therefore of 
S. If in the formula of • Article 37 we put T 3 for t and the chosen value for r, the formula 
becomes an equation giving $. The chosen value of r determines the corresponding 
values of e and § A , and the deduced value of S determines the corresponding value of s. 
Then, simultaneous values of r and s being known, all the quantities can be determined. 
It seems to be most appropriate to assume a series of suitable values of r and calculate 
the corresponding values of s. The process of finding 8, by trial, may be simplified by 
means of a theorem to the effect that the loci, in the plane of (r , s), which answer to 
constant values of t and x 09 are equally inclined to the axis of r. To prove this we 
have 

/ du\ ___ ___ j)£ Idt __ _ dXo IdXo __ (d(r\ 

\d<r/t = G<mat. dcr/ du du/ do" \ du/ XQ = const/ 

or 

' dr—ds \ ___ f dr -f- ds \ 

,dr + ds/t~ const. \dr — ds/x Q = const. 
or 

t) + (xi- = °- 

, QjT j t = const. \€vT/xq~=: const. 

This theorem shows that a point of given r on the locus * — T 3 is not far from the image 
in r = K 3 of the tangent at (R 3 , S 3 ) to the locus along which x = 0. Hence a first 
approximation to the s answering to a given r is 2S 3 -s A , where s A depends upon r in 
the known way, and therefore a first approximation to the required value of § is 

24-A. 

The junction of the second reflected wave from the left and the second middle wave 
is characterized by the value E 2 of r. If, then, the process indicated above is carried 
out for the value R 2 of r, the result is to give a pair of simultaneous values of r and s, 
which can occur in the second middle wave at the time when t = T 3 . Another pair of 
simultaneous values can be found by finding the common value of r and s which occurs 
at the central section at the same time. This is to be done by putting r = s and t = T 3 
in the formula giving t in the second middle wave, and solving the resulting equation 
for r by trial. When this is done we shall have two pairs of simultaneous values of r 
and s which occur in the left-hand half of the central part of the tube at time T 3 , and 
they are the extreme values of r and s which can occur in that part at that time. To 
obtain other pairs, we may choose an intermediate value of r, substitute in the equation 
giving t the value T 3 of * and this value of r, and find s by trial. For a first approximation 
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we may assume that the required point (r, s) is on the straight line joining the two 
extreme points whose co-ordinates have been determined previously. 

After these preliminaries the way is prepared for the numerical computation of any 
special case. 

41. Numerical Constants. — Prof. Love's investigation was undertaken in order to 
throw light on a vexed question of internal ballistics, namely, how the mass of the 
propellant should be taken into account in calculating the velocity and pressure in a 
gun. Its completion has been delayed not only by the analytical complexity of the 
problem, but also by the time required for the numerical computations. In his original 
paper Lagrange set out from a certain state of the gas assumed as a first approximation, 
namely, one in which the velocity, at a given epoch, changed uniformly from one end 
of the gas to the other. Restricting attention to the case of a very heavy gun, the 
total momentum of gas and projectile is then (M+-|C)V and the total kinetic energy 
i (M+^C)V 2 , where V is the velocity of the projectile, M its mass, and C that of the 
propelling charge. Lagrange recognized that this state of motion is dynamically 
possible only in the limiting case of small charges, but made no real progress towards 
the theory for finite charges, the development of analysis being then inadequate to the 
problem. Since the ratio C/M in modern guns, though less than with gunpowder, is 
still of the order J, the importance of a full numerical discussion of Lagrange's problem 
is evident. The calculations which follow were begun by Prof. Love, who determined 
all the fixed coefficients and the position and velocity of the projectile at various epochs. 
After verifying these figures I undertook the calculation of the distribution of pressure 
in the gas, at the times when a new type of wave was either being generated or extin- 
guished, and at the half intervals. Instantaneous combusion is assumed, as it appears, 
hopeless to attempt to allow for the gradual burning of the propellant which occurs 
in actual guns. 

It is assumed that the propellant is cordite M.D., for which the maximum pressure 
for different densities of the gas, after explosion in a closed vessel, has been measured 
by Noble.* The results at medium pressure are represented approximately by the 
formula 

#o(--l) = 9500, 

giving the pressure p in kilograms per square centimetre when Po is in absolute measure. 
This is the formula used in calculating initial pressures. The subsequent expansion of 
the gas is adiabatic, and will be represented by an equation of the form 

/I \* 
pi 1 j = const. 

\p I 
* Sir A. Noble, 'Phil. Trans./ A, vol. 205, p. 201, 1906. 
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It appears probable, for various reasons, that the mean adiabatic index y is in the 
neighbourhood of 1-2. As we are restricted to a special set of values the value 
11/9 = 1-22 is selected. 

The problem discussed in detail is that of a gun of 15 cm. calibre, mass of projectile 
50 kg. ? charge of propellant 12 kg., distance travelled by the projectile from its initial 
position of rest to the muzzle 6 metres, initial volume of gas behind the projectile 
(chamber capacity) 30 litres. It is not, of course, possible with instantaneous combus- 
tion to keep the maximum pressure the same as it would be in a gun, though the muzzle 
velocity is much the same. The maximum pressure in this case is 6333 kg. /cm. 2 . Had 
the pressure been kept down to 3000 kg. /cm. 2 by taking a smaller charge, the problem 
would have been less representative as regards muzzle velocity, and as regards the ratio 
of the masses of propellant and projectile. 

In order to exhibit both the pressure in the gun and the degree in which the back 
particles partake of the motion of the projectile, eleven planes are taken at equal distances 
apart in the undisturbed gas, the end planes coinciding with the breech and the base 
of the projectile respectively. The horizontal line at the top of Plate 1 shows their 
initial positions. These eleven planes of particles are traced throughout their motion. 
The particles originally half-way between the breech and the base of the projectile may 
be called the middle particles,* and we shall choose, as epochs for the curves of pressure 
(Plate 1), the times at which a " junction " is either at the breech or the base of the 
projectile, or has just reached the middle particles. A junction is marked with a black 
circle on the figure. 

42. Details of the Calculation (Plate 1, curve 1) (Article 10).— The early stages of 
the calculation call for no comment. We have <r = 960,536-7 cm. /sec, p = 0-4, 
p = 9500^/(1 —po) = 6333 '3 kg. /cm. 2 , c = 339*5305 cm. (the initial distance from 
the breech to the base of the projectile is \c = 169*76525 cm.), a = 177,877-1 cm. /sec, 
h = 778-0909 cm. The progressive wave which starts out from the base of the projec- 
tile reaches the middle particles (x = Jc) at time t = 0-0004772 sec. Particles between 
there and the breech are still at rest : from these particles to the base of the projectile 
the velocity of the gas increases almost uniformly to the value 99-6 m. /sec, and the 
pressure falls to 5651-3 kg. /cm. 2 . The fall of pressure is remarkable considering that 
the projectile has only moved a distance of 2-4 cm. from its initial position ; and we 
observe a finite discontinuity in the pressure gradient on the two sides of the junction. 

(Plate 1, curve 2). — The progressive wave reaches the breech at time t = 0-0009544 
sec, when the projectile has moved a distance of 9-28 cm. from its seat and has a 
velocity of 187-7 m./sec The pressure falls from 6333-3 kg. /cm. 2 at the breech to 
5097 • 2 kg. /cm. 2 at the base of the projectile. 

(Plate 1, curve 3) (Articles 12, 16-17).— The first middle wave begins at the epoch 
just mentioned, by reflexion of the progressive wave at the breech. To find when it 

* These particles must be distinguished from those of the " middle section " of the theory which here 
correspond . to the breech of the gun. 

2 H 2 
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reaches the middle particles, i.e., when the progressive wave has receded to x = Jc, 
we solve the equation (x +h) (at+h) = (h + ^c) 2 , giving £ = 0-0014785 sec. The 
velocity of the projectile at this time is 275-4 m./sec., its displacement 21-4 cm. The 
pressure falls from 5151-6 kg. /cm. 2 at the junction to 4598-7 kg./cm. 2 behind the projec- 
tile. In the first middle wave trial and error begins. At the breech u = and o- /cr 
is found by trial to give the correct value of t. For intermediate points we have theoreti- 
cally to find both u and a- by trial to make x and t correct. Actually the smallness of 
u allows us to neglect powers of uja- above the second, so that the pressure follows an 
approximately parabolic law. The difference of pressure in the first middle wave is 
quite small. At the breech we have 5170-9 kg./cm. 2 , an increase of only 19-3 kg./cm. 2 
over that at the junction, as against a drop of 552*9 kg./cm. 2 from the junction to the 
projectile. 

(Plate 1, curve 4). — The first middle wave reaches the projectile at time 
t = T x = 0-0021170 sec, when the displacement of the projectile is — X x = 42*191 cm. 
and its velocity —V x = 37175-64 cm./sec. = 371*8 m./sec. The remaining constants 
at this epoch are R x = 443,092-1, S x = 480,268-8, t x == 923,360-9. The pressure falls 
slightly from 4169-1 kg./cm. 2 at the breech to 4102-5 kg./cm. 2 at the base of the 
projectile. 

(Plate 1, curve 5) (Articles 18-25). — The first reflected wave begins at t = T t . For 
the constants we find 

log (Ao/^ 10 ) - 6*99416, log (A^ 9 ) = 9*98722, log (-2! A 2 /^ 8 ) = 7*66552, 

log (3 ! Aa/Sx 7 ) - 5-20603, log (-4 ! A^/tf) - 4*53510, log (5 ! Ag/^) - 3*35986, 

log (6 ! Ae/V) - 2*69991, log (7 ! A 7 /2i 8 ) = T'49726, log (8 ! A 8 /^ 2 ) - 0*02177, 

log (9 ! Ag/Si) - 0*29583, 

log B x = 0-33341, log B 2 = 0*84347, log B 3 = 1*66011, log B 4 - 2*49429, 

log B 5 = 3*33303, log B 6 = 418096. 

Ki = -670*58, log (~Lx) = 3*64091. 

To find when the first reflected wave reaches the middle particles, we know that 
r = E x along a junction with the first middle wave, and s is found by trial, from the 
formulae of the first middle wave, to give x = Jc. Knowing r and s, t is known : we 
find t = 0-002898 sec. The part of the first middle wave which still remains is treated 
as before. The pressure falls from 3316-0 kg./cm. 2 at the breech to 3304*3 kg./cm. 2 
at the junction. A long process is required to find the pressure in the first reflected 
wave. Writing = (r-~'R 1 )ft 1 and S = (s— S^/Si, at the base of the projectile 
= B 1 <$+B 2 <S 2 + ...+B d <$ 3 is a known function of S. We expand the formulae of 
Article 20 to give t+hja, x and (Z — E^— L^)/^ explicitly in terms of and S, and 
try different values of § until t has its required value 0-002898. Then the pressure at 
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the base of the projectile is known, and also its velocity and final position. For other 
points cp and § have to be found to make both x and t correct : the adjustment is facili- 
tated by the fact that uniform division of x corresponds approximately to uniform 
division of </:>. At the junction (x = \c) <j> = 0, and at the base of the projectile (x = 0) 
== —0-019578. Taking four values of c/> equally spaced between these, and finding 
S to give the correct t, we have four points which correspond nearly to 10, 20, 30 and 
40 per cent, division of the initial gas, and are easily adjusted to exact value by inter- 
polation. The pressure falls from 3304-3 kg. /cm. 2 at the junction to 2970*3 kg. /cm. 2 
at the base of the projectile. The projectile is displaced 75*4 cm. from its seat, and 
has velocity 466*2 m./sec. 

(Plate 1, curve 6). — The first reflected wave reaches the breech at time t = 0*003859 
sec, where the pressure is 2610*5 kg. /cm. 2 Other points are found as in the last para- 
graph. The pressure at the base of the projectile is 2161 *6 kg./cm. 2 , the displacement 
of the projectile 124*3 cm., and its velocity 550*4 m./sec. 

(Plate 1, curve 7) (Articles 30-31). — The second middle wave begins at the above 
epoch t = 0*003859, pushing back the first reflected wave along a junction s = R 1# 
This junction reaches the middle particles at time t = 0*005154. In the part of the 
first reflected wave that still remains the pressure falls from 1708*2 kg./cm. 2 at the 
junction to 1535*2 kg./cm. 2 behind the projectile. The displacement of the projectile 
is 202*1 cm. and its velocity 632*5 m./sec. The second middle wave differs from the 
first reflected wave by the presence of four additional terms with coefficients given by 

log {- % (S 2 /2E 1 )°} - 3*04397, log { m (^Rj) 7 } = 3*52835, 
log {-* 8 (S^EO 8 } = 3*60452, log { m (S 2 / 2 Ri) 9 } = 3*34841, 

where t 2 = 814,358*3 cm./sec. ; also k t = 466*85. At the breech u = or 
i p~~ $ = (S x — R^/Si, leaving <p to be found by trial. The pressure at the breech is 
1728*0 kg./cm. 2 . For intermediate points we take a number of values of <p, find S by 
trial to give the correct t, then calculate x and interpolate. 

(Plate 1, curve 8). — The second middle wave reaches the base of the projectile where 
« = S 2 = Ex, r = R 3 = 371,266*2, giving o- = S 2 , t = T 2 = 0-007137 sec. This point 
is found without trial. The pressure at the base of the projectile is 1030-2 kg./cm. 2 , 
its displacement — X 2 = 335*6 cm., and its velocity — U 2 = 71,827 cm./sec. == 
718*3 m./sec. The value of Z is Z 2 — —177*0. At the breech we have a pressure of 
1085*7 kg./cm. 2 . Other points are calculated as in the last paragraph. We have 

M^k) = 0*00034563, ?V!!i^a) = 0*0000015953, -^j^ 2 ' = - 0*000020614, 

2 2 S 2 ^2 

Fi!!i^i) = 0-00017489, Il1(|?s) = - 0*0005648, gL^|ga) = _ *004338. 
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(Plate 1, curve 9) (Articles 32-40). — Tlie second reflected wave begins at the base 
of the projectile at time T 2 and pushes back the second middle wave along a junction 
r = R 2 . Thus </) is known, and the value of S corresponding to an assigned x () is found 
by trial. We find that the junction reaches the middle particles at time t = 0*01023 
sec. The pressure at the breech is 650 - kg ./cm. 2 , at the junction 641 • kg. /cm. 2 . For 
the constants of the second reflected wave we have 



log B'x = 0*18668, log B' a = 0*27390, 
log B' B - 1*99012, log (-ci) =s 1*08789, 
logc, = 1*81519, log(-c 6 ) - 2*56985, 



logB's = 0*82262, log B' 4 =1-40104, 
logc 2 = 0*10720, log(~c 3 ) = 1*00167, 
log£ = 3*94418, log(-£ 2 ) = 5*07447, 



log ft = 6*01358, log(-&) = 6*85566, log& = 7*62931. 

The method of calculation of the pressures in the second reflected wave has been described 
in Article 40. The pressure at the base of the projectile is 581-6 kg. /cm. 2 , where the 
displacement is 571 • 9 cm. and the velocity 801 • 3 m./sec. 

The projectile is so near the muzzle at time t = 0*01023 that a fresh chart for the 
muzzle epoch (displacement 600 cm.) is unnecessary. We find for the time to the 
muzzle t = 0*01058 sec, for the muzzle velocity 807*7 m./sec, and for the pressure 
at the base of the projectile at this instant 552*6 kg. /cm. 2 . 

43. Results. — The pressure results are collected in Table I., from which Plate 1 is 
constructed. Plate 2 shows the pressures at the breech and at the base of the projec- 
tile, their ratio, the mean pressure, the displacement and velocity of the projectile, and 
a certain " energy factor " as functions of the time. The mean pressure (P in Table I.) 
is that which the cordite gases would have after adiabatic expansion, at uniform density, 
to the volume which they actually occupy at time L The work of expansion in these 
circumstances will be equal, not to the kinetic energy of the projectile, but to a greater 
kinetic energy corresponding to a fictitious mass M+aC, where 



i(M + aC)V a = fCp c 



I 



]- 



py/n- 



The " energy factor " a may be expected to vary with the distance travelled by the 
projectile : the lower values given are only approximate. 

It is difficult, after a glance at Plate 2, to resist the conclusion that the motion is 
tending to a limiting form, in which the pressure is approximately represented by 
f{y^i>{t)> with suitable functions/, <p. The energy factor a oscillates about a mean 
value of approximately 1/3, and the range of oscillation diminishes in time: similarly 
the pressure ratio oscillates about a value of approximately 0-9. Moreover, the latter 
value, like the former, can be obtained from Lagrange's approximation by suitable 
treatment.* If <p f is the pressure at the breech and p that at the base of the projectile, 



* F. Gossot and R. Liouville, * Memorial des Poudres et Salpetres,' vol. 13, p. 51, 1905. 
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p'fp is the ratio of the momenta of gun and projectile, that is (M+-|C)/M, so that the 
pressure ratio is approximately — = 0-893. The agreement is to be expected ; 

for Table I. shows how little, relatively, yfy varies with y 09 so that Lagrange's approxi- 
mation leads to little error in the total energy and momentum. 

44. Calculation of Recoil. — Prof. Love's theory also enables us to calculate the distance 
recoiled by a very heavy gun while the projectile is travelling to the muzzle : this is 
important since the distance can also be found experimentally. We take from Plate 2 

the values of p ! and p at intervals of 0-0005 sec. to the muzzle, and calculate p' dt 

and \pdt by approximate integration. These quantities are proportional to M'V 

and MV, where M 7 , M are the masses of gun and projectile and V, V their velocities. 
A second integration gives M'S' and MS, where S' and S are the distances travelled by 
gun and projectile. For the muzzle epoch we find, in the present problem, M'S'/MS = 
1/0-879. The recoil distance S' of the gun is therefore the same as for a massless pro- 
pellant and a projectile of mass M/0- 879 = 56-9 kg., an addition of 0-57 times the 
mass of the propellant to that of the projectile. Lagrange's approximation gives 
• 5. Crantz* measured the recoil distance of a rifle, with comparatively slow combustion 
of the propellant, and obtained factors 0-496, 0-497, 0-477, mean 0-493. The theory 
of limiting motion would seem to apply with almost equal force to the case of slow 
combustion ; and thus we may regard Cranz's experiment as confirming the recoil 
factor | and therefore (indirectly) the energy factor 1/3. Prof. Love has worked out 
the energy factor for a light projectile of mass 25 kg., and 12 kg. propellant, at epochs 
corresponding to (4) and (8) in Table I. The values are 0-335 and 0-333. 

45. A Special Solution of the Hydrodynamical Equations. — Prof. Love's theory having 
suggested the possibility of the motion tending to a limiting form, it remains to show 
that the hydrodynamical equations admit of a particular solution in which the pressure 
is of the form f(y ) <p (t). We shall see that the pressure ratio and energy factor corre- 
sponding to this exact solution agree closely with those already calculated, and thus 
support is lent to the view that the limiting motion would be developed sooner or later 
with other initial conditions, e.g., with gradual introduction of gas from a burning 
propellant. If y is, as above, the initial distance of a particle from the breech and y 
its distance at time t, the general hydrodynamical equation is 



3V /I A y (1 dy A-y- 1 1 3V 

ot \p I \p dy J p dy 



2 



Write temporarily x = y , z = y—p y . Then 

df Po [ P0} dx 2 \ dx, 
* C. Ceanz, ' Zeitschr. f. d. ges. Schiess-u. Sprengstoffwesen/ vol. 2, p. 345, 1907, 
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A solution of the type z — f(x) <j> (t) is possible if and only if 

/' (x) = a/(x) {/' (x)}y-\ {<p (t)}y +» (t) = B, 

where A and B are constants connected by the equation 

B = ?& {1 - poY A. 

Po 

If S is the area of the cross-section, the equation of motion of the projectile, which is 
supposed to be at x = b, is p$ = M/(6) $" (t). Now in general 

p=po(i-poy{f(x)<p(t)}~\ 

Hence the equation of motion of the projectile is satisfied if 

MB/(&) = Sp (l-po)M/(&)}- 7 ? 



or 

A = 



ybf{b){f(b)Y' 



where € = C/M = S&p /M is the ratio of the mass of the propellant to that of the projec- 
tile. Writing w =f(x) and q = dw/dx, the first integral of the differential equation 
for w is 

7-1 - 2 1 

q ~~ A( r -l)a 2 -w 2 ' 

where a is a constant. Since f(x) vanishes with x 9 the final integral is 

i 

% io ._!__ f o 1 y— l 







(a 2 —w 2 )y- 1 chv = <^ -r-7 7? f #. 



Writing c =f(b) for the length of the column of gas at the instant considered, we have 

therefore 

1 

2 |7-i 7 r /v/tMv-i 2 1 



Substituting for A we have 

' (6] " 2y&C A(y-l) I 2y6c J ^ ° j ' 

so that • 

V y 2yC 

This equation determines c/a, and when it is known w is given by 



% w 



x 2yC O 
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The pressure ratio between the two ends of the gas is E = {f'(0)/f'(b)} y , where 
if {VlfiO)} 1 ' 1 = <*'/(<*'-<?). Hence 

Y 
.2 ^ 



Jti = 



a 2 



Writing c = a sin 6 we find R = cos m 0, where m = 2y/(y — l) and is found by trial 

from the equation 

sin "* 



.w 



cos" 1 " 1 0cZ0 = 



cos J o m 

In the case of y = 11/9, m = 11, we find, after some analytical reduction, the 
expansion 

"R — 1 — -L_L 7 £2_ 14 913 3 , 

X\, __ i— 2 e-h^ 64 e — 8TT20 6 + ••• j 

valid for small values of <?. Either method gives R = 0-894 for e = 12/50, the corre- 
sponding value of being 8° 9' 6". It will be noticed that although R is not equal to 
(l+i e ) _1 to the second order, the approximation is still a remarkably good one. The 
present theory will appear more satisfactory, as it is based on an exact solution valid 
for all values of e. 

As regards the energy factor, the previous definition in terms of the work done from 
an initial state of uniform density is not convenient, as this state is not one of the previous 
states of the gas. We may, however, define the energy factor in such a way that the 
kinetic energy of the gas is ae times that of the projectile.* 

Corresponding to the initial distance x from the breech we have in general w = a sin 0, 

where 

§ _ x 

cos m X (hdd> = K r , 


and 

K = f^os™" 1 0cZ0 = ^p!^. 

Jo m sin 

The corresponding velocity is V sin </>/sin 0. If x-{-dx corresponds to (p J r d<p, 
K.dxfb = cos m-1 <pd(p. The kinetic energy of the gas is 



x CV 2 rd 
2 K sin 2 



cos™ 1 <p sin 2 (j>d<j), 



and that of the projectile |-MV 2 . Hence by definition 



m 

q ■ == __ — 

e cos™ sin 



cos m 1 (p sin 2 cj>d(j>, 



* This was not done above because the problem would naturally present itself in the other form in 
practical calculations, where we should seek a factor which will make the kinetic energy of the projectile 
equal to the work of an assumed massless propellant. 
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where 8 is given by the equation already written down. Using the reduction 

formulae we find 

e— m sin 2 

(m+l) esin 2 #' 

giving a = 0*325 when m = 11 and e = 12/50. The expansion formula, to the first 
power in e, is 

a = £— rib <?• 

46. Application to Ballistics. — To resume, Prof. Love's theory supports the factors 
-| and |- up to considerable values of C/M, and shows further that the ratio of the pressures 
on projectile and breech (Plate 2) begins at once to oscillate about its mean value, 
reaching its first minimum when the projectile has travelled a distance of only two-thirds 
of a calibre. We may remark that no support is lent to the theory which appears to be 
favoured by Charbonnier* of more or less violent impulses of pressure on the base 
of the projectile : the discontinuity is at most one of pressure gradient, which becomes 
less and less as the motion proceeds. What would happen with gradual introduction 
of gas from a burning propellant is more conjectural, but nevertheless it seems of interest 
to examine the consequences of the assumption that the limiting state of motion, contem- 
plated above, is developed almost at once, and maintained ever after. The considerations 
which we shall advance have no pretence to rigour, f 

It is usual to measure maximum pressures in guns by crusher gauges placed at or near 

the breech. Let P be the pressure at the breech, P(l— C/2M) that at the base of the 

projectile at any time, powers of C/M above the first being neglected. Compare the 

actual motion with that for a massless gas of the same thermodynamical properties, 

and a projectile of mass m. Then for identical motion of the two projectiles, with 

m/M - 1+C/3M, 

p __ m 

P(1-C/2M) ~~M ? 

or p/T = 1— -C/6M. In order to keep up the parallelism of motion we have to ensure 

that equal quantities of propellant are burnt in equal times. The rate of regression of 

the surface ol colloidal propellants at different pressures has been measured by Vieille 

in a famous research.*! Mansell, who examined cordite M.D. by Vieille ? s method,§ 

found a rate of regression in a closed vessel approximately proportional to the pressure. 

If D and d are the diameters of cordite in the two cases (or more generally numbers 

proportional to the linear dimensions of the grain), equal generation of gas corresponds 

approximately to the condition 

d_ __ p 

D~ P(1-C/4M) 5 

* P. Charbonnier, ' Traite de Balistique Interieure,' Paris, 0. Doin, p. 91. 

t See also F. Gossot and E. Liouville, loc. cit, pp. 50-58 ; vol. 17, pp. 61-66, 1914. 

% P. Vieille, ' Memorial des Poudres et Salpetres,' vol. 6, p. 256, 1893. 

§ J. H. Mansell, ' Phil. Trans.,' A, vol. 207, p. 243, 1908. 
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since P(l — -C/4M) is the mean pressure in the. first gun. Hence we find that the relation 
of the ballistic constants in the two guns is 

w 1 , C d 1 , P , , C 



W 3M 5 D 12M ? p ' 6M' 

The first two equations give the projectile and the size of cordite to be used in the ideal 
calculation, the third equation the ratio in which the calculated maximum pressure is 
to be increased. The theory would be seriously invalidated if the ratio of the pressures 
recorded by crusher gauges in the breech and in the base of the projectile is not approxi- 
mately 1+C/2M, and does not apply to extraordinary experiments with very quick 
combustion (for which the pressure ratio is, of course, nearly 1). 

From a few calculations I have made with full charges in guns, it appears that the 
empirical rule of adding one-half of the mass of the propellant to the mass of the projec- 
tile, without other change of ballistic constant, gives approximately correct pressures, 
while muzzle velocities are about If per cent. low. The fraction one-third gives approxi- 
mately correct muzzle velocities, but maximum pressures about 4 per cent. low. 

An important factor in the future progress of internal ballistics would seem to be 
the determination of the rate of regression of colloidal propellants as a function of both 
temperature and pressure. Hitherto only the latter has been taken into account, 
although some experiments of Wolff* show a falling off of the rate for small charges 
in a very small closed vessel, which appears to be due to loss of temperature. The effect 
would be enhanced in a gun, where the whole mass of gas is cooled by expansion, instead 
of being cooled relatively strongly near the surface. One consequence of diminished 
burning would be the occurrence of unconsumed cordite at velocities higher than those 
which Sebert and Hugoniot's formula would give with a burning constant derived 
from experiments in closed vessels. 



* W. Wolff, ' Kriegstechniselie Zeitschr.,' vol. 6, p. 1, 1903. 
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